SMOOTH DYNAMICS AND 
NEW THEORETICAL IDEAS IN NONEQUILIBRIUM STATISTICAL MECHANICS. 



by David Ruelle* 



Abstract. This paper reviews various applications of the theory of smooth dynamical 
systems to conceptual problems of nonequilibrium statistical mecanics. We adopt a new 
point of view which has emerged progressively in recent years, and which takes seriously 
into account the chaotic character of the microscopic time evolution. The emphasis is on 
nonequilibrium steady states rather than the traditional approach to equilibrium point of 
view of Boltzmann. The nonequilibrium steady states, in presence of a Gaussian thermo- 
stat, are described by SRB measures. In terms of these one can prove the Gallavotti-Cohen 
fluctuation theorem. One can also prove a general linear response formula and study its 
consequences, which are not restricted to near equilibrium situations. At equilibrium one 
recovers in particular the Onsager reciprocity relations. Under suitable conditions the 
nonequilibrium steady states satisfy the pairing theorem of Dettmann and Morriss. The 
results just mentioned hold so far only for classical systems; they do not involve large size, 
i.e. they hold without a thermodynamic limit. 
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Chapter I. A GENERAL DISCUSSION. 



1. Introduction. 

Statistical mechanics, which was created at the end of the 19-th century by such peo- 
ple as Maxwell, Boltzmann, and Gibbs, consists of two rather different parts: equilibrium 
and nonequilibrium statistical mechanics. The success of equilibrium statistical mechanics 
has been spectacular. It has been developed to a high degree of mathematical sophis- 
tication, and applied with success to subtle physical problems like the study of critical 
phenomena. Equilibrium statistical mechanics also has highly nontrivial connections with 
the mathematical theory of smooth dynamical systems and the physical theory of quantum 
fields. 

By contrast, the progress of nonequilibrium statistical mechanics has been much 
slower. We still depend on the insights of Boltzmann for our basic understanding of 
irreversibility, and this understanding remains rather qualitative. Further progress has 
been mostly on dissipative phenomena close to equilibrium: Onsager reciprocity rela- 
tions, Green-Kubo formula, and related results. Yet, there is currently a strong revival of 
nonequilibrium statistical mechanics, based on taking seriously the complications of the 
underlying microscopic dynamics. The purpose of the present review paper is to discuss 
some of the new ideas in the dynamical approach to the problems of nonequilibrium. For 
a different presentation, see [28] . 

This paper has grown out of lectures given in the Mathematics Department at Rutgers 
University in 1997 and 1998, where I benefitted from a very interactive and exceptionally 
competent audience. I have thus used in what follows what I learnt from discussions with 
Joel Lebowitz, Giovanni Gallavotti, Sheldon Goldstein, Eugene Speer, Michael Kiesling, 
Frederico Bonetto, and others. Let all of them be thanked. 

2. Irreversibility. 

The rest of this chapter is devoted to a discussion of irreversibility based on Boltz- 
mann's ideas. This will serve as an introduction and background for our later discussions. 

It is a fact of common observation that the behavior of bulk matter is often irreversible: 
if a warm and a cold body are put in contact, they will equalize their temperatures, but 
once they are at the same temperature they will not spontaneously go back to the warm 
and cold situation. Such facts have been organized into a body of knowledge named ther- 
modynamics. According to thermodynamics, a number called entropy can be associated 
with macroscopic systems which are (roughly speaking) locally in equilibrium. The defi- 
nition is such that, when the system is isolated, its entropy can only increase with time or 
stay constant. (This is the celebrated second law of thermodynamics) . A strict increases in 
entropy corresponds to an irreversible process. Such processes are also called dissipative, 
because they dissipate noble forms of energy (like mechanical energy) into heat. The flow 
of a viscous fluid, the passage of an electric current through a conductor, or a chemical 
reaction like a combustion are typical dissipative phenomena. The purpose of nonequilib- 
rium statistical mechanics is to explain irreversibility on the basis of microscopic dynamics, 
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(2) 



and to give quantitative predictions for dissipative phenomena. In what follows we shall 
assume that the microscopic dynamics is classical. Unfortunately we shall have little to 
say about the quantum case. 

3. Dynamics and entropy. 

We shall begin with a naive discussion, and then see how we have to modify it to 
avoid the difficulties that will arise. The microscopic time evolution (/*) which we want 
to consider is determined by an evolution equation 

(D 

in a phase space M. More precisely, for an isolated system with Hamiltonian H, we may 
rewrite (1) as 

±(p\ = /-0 q tf(p,q 
dt\q) \d p H(p,q) 

where p, and q are iV-dimensional for a system with N degrees of freedom ( Ai is thus 2N- 
dimensional) . Note that we are interested in a macroscopic description of a macroscopic 
system where N is large, and microscopic details cannot be observed. In fact many different 
points (p, q) in phase space have the same macroscopic appearance. It appears therefore 
reasonable to describe a macroscopic state of our system by a probability measure m(dx) 
on A4. At equilibrium, m(dx) should be an (/^-invariant measure. 

Remember now that the Hamiltonian time evolution determined by (2) preserves the 
energy and the volume element dpdq in phase space. This leads to the choice of the 

microcanonical ensemble 

m K {dpdq) = -^-<y(fl"(p,q) - K)dpdq (3) 

to describe the equilibrium state of energy K. [Qk is a normalization constant, and the 
ergodic hypothesis asserts that tuk is ergodic]. Note that the support of tuk is the energy 
shell 

A^ = {(p,q):tf(p,q)=K} 

If the probability measure m(dpdq) has density m(p, q) with respect to dpdq, we 
associate with the state described by this measure the entropy 

S(m) = - J dpdqm(p,q)logm(p,q) (4) 

This is what Lebowitz [62] calls the Gibbs entropy, it is the accepted expression for the 
entropy of equilibrium states*, and it remains reasonable outside of equilibrium, as infor- 
mation theory for instance would indicate. Writing x instead of (p, q), we know that the 
density of m = f^m is given by 

m{.f~ tx ) 



m(x) = 



* There is a minor technical complication here. While (4) gives the right result for 
the canonical ensemble one should, for the microcanonical ensemble, replace the reference 
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where J t is the Jacobian determinant of /*. We have thus 

S(m) = — J dx m(x) log m(x) (5) 

S(m) = - / dxm(x) \ogm(x) = - dx . , , log — , , 

dxm(x) log ~~ = S'(m) + / dxm(x) log Jt(x) (6) 
•'/(•<•) ./ 

In the Hamiltonian situation that we are considering for the moment, since the volume 
element is preserved by /*, J t = 1, and therefore S(m) = S(m). 

So, we have a problem: the entropy seems to remain constant in time. In fact, we 
could have expected that there would be a problem because the Hamiltonian evolution (2) 
has time reversal invariance (to be discussed later), while we want to prove an increase of 
entropy which does not respect time reversal. We shall now present Boltzmann's way out 
of this difficulty. 

4. The classical ideas of Boltzmann and his followers. 

Let us cut M. into cells Cj so that the coordinates p, q have roughly fixed values 
in a cell. In particular all points in Cj are macroscopically equivalent (but different Cj 
may be macroscopically indistinguishable). Instead of describing a macroscopic state by 
a probability measure m, we may thus give the weights m(c^). Now, time evolution will 
usually distort the cells, so that each /*Cj will now intersect a large number of cells. If the 
initial state m occupies M cells with weights 1/M (taken to be equal for simplicity), the 
state /**m will occupy (thinly) Af f cells with weights 1/A/"*, where N l may be much larger 
than M. If the Cj have side h, we have 

log J\f =— ^ m(ci) log m(ci) — / dx m(x) log(m(x). h 2N ) = S(m) — 2N log h 
i J 

i.e., the entropy associated with m is roughly log A/" + 2Nlogh. The apparent entropy 
associated with /**m is similarly logA/"* + 2iVlog/i; it differs from S(m) because the 
density m, which may fluctuate rapidly in a cell, has been replaced by an average for the 
computation of the apparent entropy. The entropy increase logA/"* — log M is due to the 
fact that the initial state m is concentrated in a small region of phase space, and becomes 
by time evolution spread (thinly) over a much larger region. The time evolved state, after a 



measure dpdq in (3) by S(H(p, q) — K)dpdc[. This point is best appreciated in the 
light of the theory of equivalence of ensembles in equilibrium statistical mechanics. For 
our purposes, the easiest is to replace 5 in (3) by the characteristic function of [0, e] for 
small e. The measure uik is then still invariant, but no longer ergodic. Note also that the 
traditional definition of the entropy has a factor k (Bolzmann's constant) in the right-hand 
side of (4). 
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little smoothing {coarse graining), has a strictly larger entropy than the initial state. This 
gives a microscopic interpretation of the second law of thermodynamics. In specific physical 
examples (like that of two bodies in contact with initially different temperatures) one sees 
that the time evolved state has correlations (say between the microscopic states of the two 
bodies, after their temperatures have equalized) which are macroscopically unobservable. 
In the case of a macroscopic system locally close to equilibrium (small regions of space 
have a definite temperature, pressure,. . . ) the above classical ideas of Boltzmann have been 
expressed in particularly clear and compelling manner by J. L. Lebowitz [62]. He defines a 
local Boltzmann entropy to be the equilibrium entropy corresponding to the temperature, 
pressure,. .. which approximately describe locally a nonequilibrium state. The integral 
over space of the Boltzmann entropy is then what we have called apparent entropy, and 
is different from the Gibbs entropy defined by (4). These ideas would deserve a fuller 
discussion, but here we shall be content to refer to [62], and to Boltzmann's original works. 

There is still some opposition to Boltzmann's ideas (notably by I. Prigogine and his 
school), but most workers in the area accept them, and so shall we. We shall however 
develop a formalism which is both rather different from and completely compatible with 
the ideas of Boltzmann and Lebowitz just discussed. 

5. The description of states by probability measures. 

In the above discussion, we have chosen to describe states by probability measures. 
One may object that the state of a (classical) system is represented by a point in phase 
space rather than by a probability measure. But for a many-particle system like those 
of interest in statistical mechanics it is practically impossible to know the position of all 
the particles, which changes rapidly with time anyway. The information that we have 
is macroscopic or statistical* and it is convenient to take as description of our system a 
probability distribution p on phase space compatible with the information available to us. 
Trying to define this p more precisely at this stage leads to the usual kind of difficulties 
that arise when one wants to obtain from a definition what should really come as a result 
of theory. 

For our later purposes it is technically important to work with states described by 
probability measures. Eventually, we shall obtain results about individual points of phase 
space (true almost everywhere with respect to some measure). From a physical point of 
view, it would be desirable to make use here of points of phase space which are typical for 
a macroscopic state of our system. Unfortunately, a serious discussion of this point seems 
out of reach in the present state of the theory. 

6. Beyond Boltzmann. 

There are two reasons why one would want to go beyond the ideas presented above. 
One concerns explicit calculations, like that of a rate of entropy production; the other 
concerns situations far from equilibrium. We discuss these two points successively. 

If one follows Boltzmann and uses a decomposition of phase space into cells to compute 
entropy changes, the result need not be monotone in time, and will in general depend on 

* Information on local microscopic correlations is obtained by diffraction techniques. 
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the particular decomposition used. Only after taking a limit t — > oo can one let the size 
of cells tend to 0, and obtain a result independent of the choice of coarse graining. This 
leaves open the problem of computing the entropy production per unit time. 

The idea of using local Boltzmann entropies works only for a macroscopic system 
locally close to equilibrium, and one may wonder what happens far from equilibrium. In 
fact one finds statements in the literature that biological processes are far from equilibrium 
(which is true), and that they may violate the second law of thermodynamics (which is not 
true). To see that life processes or other processes far from equilibrium cannot violate the 
second law, we can imagine a power plant fueled by sea water: it would produce electric 
current and reject colder water in the sea. Inside the plant there would be some life form or 
other physico-chemical system functioning far from equilibrium and violating the second 
law of thermodynamics. The evidence is that this is impossible, even though we cannot 
follow everything that happens inside the plant in terms of Boltzmann entropies of systems 
close to equilibrium. In fact, Boltzmann's explanation of irreversibility reproduced above 
applies also here, and the only unsatisfactory feature is that we do not have an effective 
definition of entropy far from equilibrium. 

The new formalism which we shall introduce below is quite in agreement with the ideas 
of Boltzmann which we have described. We shall however define the physical entropy by 
(4) (Gibbs entropy). To avoid the conclusion that this entropy does not change in time 
for a Hamiltonian time evolution, we shall idealize our physical setup differently, and in 
particular introduce a thermostat. 

7. Thermostats. 

Instead of investigating the approach to equilibrium as we have done above following 
Boltzmann, we can try to produce and study nonequilibrium steady states. To keep a finite 
system outside of equilibrium we subject it to non-Hamiltonian forces. We consider thus 
an evolution of the form (1), but not (2). Since we no longer have conservation of energy, 
x(t) cannot be expected to stay in a bounded region of phase space. This means that the 
system will heat up. Indeed, this is what is observed experimentally: dissipative systems 
produce heat. An experimentalist will eliminate excess heat by use of a thermostat, and 
if we want to study nonequilibrium steady states we have to introduce the mathematical 
equivalent of a thermostat. 

In the lab, the system in which we are interested (called small system) would be cou- 
pled with a large system constituting the thermostat. The obvious role of the large system 
is to take up the heat produced in the small system. At the same time, the thermostat 
allows entropy to be produced in the small system by a mechanism discussed above: mi- 
croscopic correlations which exist between the small and the large system are rendered 
unobservable by the time evolution. An exact study of the pair small+large system would 
lead to the same problems that we have met above with Boltzmann's approach*. Note 

* For such an approach, see Jaksic and Pillet [53] , where the large system has infinitely 
many noninteracting degrees of freedom. Studying the diffusion and loss of correlations 
in an infinite interacting system (say a system of particles with Hamiltonian interactions) 
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however that the thermostats used in the lab are such that their state changes as little 
as possible under the influence of the small system. For instance the small system will 
consist of a small amount of fluid, surrounded by a big chunk of copper constituting the 
large system: because of the high thermal conductivity of copper, and the bulk of the 
large system, the temperature at the fluid-copper interface will remain constant to a good 
precision. In conclusion, the experimentalist tries to build an ideal thermostat, and we 
might as well do the same. We replace thus (1) by 

(IT 

- = F(x) + Q(u(t),x) 

where the effect Q(u(t), x) of the thermostat depends on the state x of the small system, 
and on the state u(t) of the thermostat, but the time evolution t — > u(t) does not depend 
on the state x of the small system. 

We may think of ui{t) as random (corresponding to a random thermostat), but the 
simplest choice is to take to constant, and use Q(x) = Q(u, x) to keep x(t) on a compact 
manifold M. For instance if M is the manifold {x : h(x) = K} : we may take 

rv \ (F(x)-gradh(x)) 
(gradn,(x) ■ granh(x)) 

This is the so-called Gaussian thermostat [51], [32]. We shall be particularly interested later 
in the isokinetic thermostat, which is the special case where x = (p, q) and h(x) = p- p/2m 
(kinetic energy). 

8. Nonequilibrium steady states. 

We assume for the moment that the phase space of our system is reduced by the action 
of a thermostat to be a compact manifold M. The time evolution equation on M has the 
same form as (1): 

I - Fix) (7) 

where the vector field F on M now describes both the effect of nonhamiltonian forces 
and of the thermostat. Note that (7) describes a general smooth evolution, and one may 
wonder if anything of physical interest is preserved at this level of generality. Perhaps 
surprisingly, the answer is yes, as we see when we ask what are the physical stationary 
states for (7). 

We start with a probability measure m on M such that m(dx) = m{x) dx, where dx is 
the volume element for some Riemann metric on M (for simplicity, we shall say that m is 



appears to be very difficult in general, because the same particles may interact again and 
again, and it is hard to keep track of the correlations resulting from earlier interactions. 
This difficulty was bypassed by Lanford when he studied the Boltzmann equation in the 
Grad limit [60], because in that limit, two particles that collide once will never see each 
other again. 
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absolutely continuous) . At time t, m becomes (/*)*m, which still is absolutely continuous. 
If (f t )*m has a limit p when t — > oo, then p is invariant under time evolution, and in 
general singular with respect to the Riemann volume element dx. [A time evolution of 
the form (7) does not in general have an absolutely continuous invariant measure]. The 
probability measures 



or more generally 



p = lim (f )*m 

t^oo 



1 '* 



lim - / dr{f T Ym (8) 

t^oo t ' 







(with m absolutely continuous) are natural candidates to describe nonequilibrium station- 
ary states, or nonequilibrium steady states. Examples of such measures are the SRB states 
(see [92], [79], [14], [66], [67]), which will be discussed in greater detail in Chapter 2. 

9. Entropy production. 

We return now to our calculation (6) of the entropy production: 



S(m) — S(m) = J dxm(x)\ogJ t (x) 

where m is the density of m = f u m. This is the amount of entropy gained by the system 
under the action of the external forces and thermostat in time t. The amount of entropy 
produced by the system and given to the external world in one unit of time is thus (with 
J = Ji) 

e/(m) = — J m(dx) log J(x) 

Notice that this expression makes sense also when m is a singular measure. The average 
entropy production in t units of time is 

\^e f {f k m) = \[S{m) - S{m)\ (9) 

k=0 

When t — > oo, this tends according to (8) towards 

e f (p) = -J p(dx) log J(x) (10) 

which is thus the entropy production in the state p (more precisely, the entropy production 
per unit time in the nonequilibrium steady state p). Using (7) we can also write 

e, W =-/^)divF W (H) 



Notice that the entropy S in (9) is bounded above (see the definition (5)), so that 
e/(p) > 0, and in many cases e/(p) > as we shall see later. Notice that for an arbi- 
trary probability measure p (invariant or not), e/(p) may be positive or negative, but the 
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definition (8) of p makes a choice of the direction of time, and results in positive entropy 
production. It may appear paradoxical that the state p, which does not change in time, 
constantly gives entropy to the outside world. The solution of the paradox is that p is 
(normally) a singular measure and therefore has entropy — oo: the nonequilibrium steady 
state p is thus a bottomless source of entropy. 

10. A new idealization of nonequilibrium processes. 

We have now reached a new framework idealizing nonequilibrium processes. Instead 
of following Boltzmann in his study of approach to equilibrium, we try to understand the 
nonequilibrium steady states as given by (8). For the definition of the entropy production 
ef(p) we use the rate of phase space contraction (10) or (11). (An early reference proposing 
this definition is Andrey [1]). In Chapter 2, we shall discuss the SRB states, which provide 
a mathematically precise definition of the nonequilibrium steady states. After that, the 
idea is to use SRB states to make interesting physical predictions, making hyperbolicity 
assumptions (this will be explained later) as strong as needed to get results. Such a 
program was advocated early by Ruelle, but only recently were interesting results actually 
obtained, the first one being the fluctuation theorem of Gallavotti and Cohen [44], [45]. 
This will be discussed in Chapter 3. In later chapters we shall review the growing body of 
knowledge which is emerging in this new approach to nonequilibrium statistical mechanics. 

It should be stressed at this point that the program just outlined is strictly classical, 
and has for the moment no quantum version. Note in this respect that the classical entropy 
may go to — oo, while the quantum entropy is defined to be > 0, this makes it difficult 
to propose a quantum nonequilibrium steady state which would be a bottomless source of 
entropy. 

11. The diversity of nonequilibrium regimes. 

Many important nonequilibrium systems are locally close to equilibrium, and the clas- 
sical nonequilibrium studies have concentrated on that case, yielding such results as the 
Onsager reciprocity relations and the Green-Kubo formula. Note however that chemical 
reactions are often far from equilibrium. More exotic nonequilibrium systems of inter- 
est are provided by metastable states. Since quantum measurements (and the associated 
collapse of wave packets) typically involve metastable states, one would like to have a rea- 
sonable fundamental understanding of those states. Another class of exotic systems are 
spin glasses, which are almost part of equilibrium statistical mechanics, but evolve slowly, 
with extremely long relaxation times. 
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Chapter 2. SRB STATES. 



1. Further discussion of nonequilibrium steady states. 

In Chapter 1 we have proposed a definition (see (1.8)) for nonequilibrium steady states 
p. We now make this definition more precise and analyze it further. Write 



1 '* 



p = w.lim^oo - dr (/ r )*m (1) 
f Jo 

m a.c. probability measure (2) 

p ergodic (3) 

In (1), w.lim is the weak or vague limit defined by (w.limm t )($) = lim(mt($)) for all 
continuous $ : M — > C. The set of probability measures on M is compact and metrizable 
for the vague topology. There are thus always sequences (£/-) tending to oo such that 

1 f tk 

p = w.hnifc^oo — / dr (f T )*m 
tk Jo 

exists; p is automatically (/ T )* invariant. By (2), we ask that the probability measure 
m be a.c. (absolutely continuous) with respect to the Riemann volume element dx (with 
respect to any metric) on M. The condition (3) is discussed below. 

Physically, we consider a system which in the distant past was in equilibrium with 
respect to a Hamiltonian time evolution 

and described by the probability measure m on the energy surface M. According to 
conventional wisdom, m is the microcanonical ensemble, which satisfies (2), and is ergodic 
with respect to the time evolution (4) when the ergodic hypothesis is satisfied*. For our 
purposes, we might also suppose that m is an ergodic component of the microcanonical 
ensemble or an integral over such ergodic components, provided (2) is satisfied. 

We assume now that, at some point in the distant past, (4) was replaced by the time 
evolution 

representing nonhamiltonian forces plus a thermostat keeping x on M; we write the general 
solution of (5) asxw f l x. We are interested in time averages of f l x for m-almost all x. 
Suppose therefore that 

1 rt 



w.lim^oo - / dr5fT X (6) 



o 



* Even if the ergodic hypothesis is accepted, the physical justification of the micro- 
canonical ensemble remains a delicate problem, which we shall not further discuss here. 
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exists for m-almost all x. In particular, with p denned by (1), we have 



p = J m {dx) p x (7) 

If (3) holds*, (7) is equivalent to 

p x = P for m-almost all x 

[<= is obvious; if =>- did not hold (7) would give a non-trivial decomposition p = ap' + (1 — 
a)p" in contradiction with ergodicity]. 

As we have just seen, the ergodic assumption (3) allows us to replace the study of (6) 
by the study of (1), with the condition (2). This has interesting consequences, as we shall 
see, but note that (1-3) are not always satisfyable simultaneously [consider for instance 
the case F = 0]. To study nonequilibrium steady states we shall modify or strengthen 
the conditions (1-3) in various ways. We may for simplicity replace the continuous time 
t G R by a discrete time t G Z. If we assume uniform hyperbolicity (Sections 2 and 3), the 
states satisfying (1-3) are called SRB states (Section 4). More general SRB states have 
been studied by Ledrappier, Strelcyn, and Young (Section 5) and satisfy the conditions 
(1) under a weak (nonuniform) hyperbolicity condition. 

2. Uniform hyperbolicity: Anosov diffeomorphisms and flows. 

Let M be a compact connected manifold. In what follows we shall be concerned with 
a time evolution (/*) which either has discrete time t G Z, and is given by the iterates of a 
diffeomorphism f of M, or has continuous time t G R, and is the flow generated by some 
vector field F on M. We assume that / or F are of class C 1+a (Holder continuous first 
derivatives). 

The Anosov property is an assumption of strong chaoticity (in physical terms) or 
uniform hyperbolicity (in mathematical terms**). Choose a Riemann metric on M. The 
diffeomorphism / is Anosov if there are a continuous T/-invariant splitting of the tangent 
bundle: TM = E s © E u and constants C > 0, 6 > 1 such that 

[V - UJ \\Tf-^\\ <ce- 1 if ieE u [6) 

One can show that x i— > E*, E% are Holder continuous, but not C 1 in general. The flow 
(/*) is Anosov if there are a continuous (Tf*)-invariant splitting*** TM = H.FQ)E S @E U 
and constants C > 0, 9 > 1 such that (8) again holds. 



* Suppose that p is not ergodic but that p = ap' + (1 — a)p" with p' ergodic and a ^ 0. 
Writing S = {x : p x = p'}, we have m(S) = a and p' = f m'(dx) p x with m! = a _1 X5.m. 
Therefore, (1-3) hold with p, m replaced by p', m' . 

** For background see for instance Smale [93] or Bowen [13] 
*** Remember that F is the vector field generating the flow: F(x) = df l x/dt\t=Q. There- 
fore R.F is the one-dimensional subbundle in the direction of the flow. 
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In what follows we shall assume that the periodic orbits are dense in M. This is 
conjectured to hold automatically for Anosov diffeomorphisms, but there is a counterex- 
ample for flows (see [34]). Since M is connected, we are thus dealing with what is called 
technically a mixing Anosov diffeomorphism /, or a transitive Anosov flow (/*). There is a 
powerful method, called symbolic dynamics, for the study of Anosov diffeomorphisms and 
flows. Symbolic dynamics (see Sinai [90]) is based on the existence of Markov partitions 
(Sinai [91], Ratner [76]). Another, more geometric approach based on shadowing (Bowen 
[13]) is also available, and will be used in Chapter 3. 

3. Uniform hyperbolicity: Axiom A diffeomorphisms and flows. 

Smale [93] has achieved an important generalization of Anosov dynamical systems by 
imposing hyperbolicity only on a subset O (the nonwandering set) of the manifold M. 
A point x G M is wandering if there is an open set O 3 x such that O H fO 7^ for 
\t\ > 1. The points of M which are not wandering constitute the nonwandering set O. A 
diffeomorphism or flow satisfies Axiom A if the following conditions hold 

(Aa) there is a continuous (T/*)-invariant splitting of TqM (the tangent bundle re- 
stricted to the nonwandering set) verifying the hyperbolicity conditions of Section 2 above. 

(Ab) the periodic orbits are dense in Q. 

Under these conditions, O is a finite union of disjoint compact (/*) -invariant sets B 
(called basic sets) on which (/*) is topologically transitive*. This result is known as Smale's 
spectral decomposition theorem. 

If there is an open set U D B such that C\t>ofU = £?, the basic set B is called an 
attractor. The set {x G M : limt^ +00 d(fx, B) = 0} = Utyof^U is the basin of attraction 
of the attractor B. 

Let B be a basic set. Given x G B and e > 0, write 



Then for sufficienty small e, VJ e and V" e are pieces of smooth manifolds, called local stable 
and local unstable manifold respectively**, and tangent to E*, respectively. 

* (/*) is topologically transitive on B if there is x G B such that the orbit (f l x) is dense 
in B. 

** There are also global stable and unstable manifolds defined by 



V s xe = {yeM: d(fy, fx) < e for t > 0, and lim t 



+00 d(f t y,f t x) = 0} 



V£ e = {yeM: d(fy, fx) < e for t < 0, and lim t 



—00 



d(fy,fx) = 0} 



V s x = {yeM: lim d(fy,fx) = 0} 
t^+00 



VI = {y G M : lim d(fy, fx) = 0} 
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A basic set B is an attractor if and only if the stable manifolds V| e for x G B cover a 
neighborhood of B. Also, a basic set i? is an attractor if and only if the unstable manifolds 
V" e for x E B are contained in B (see [14]). 

Markov partitions, symbolic dynamics (see Bowen [9], [12]), and shadowing (see Bowen 
[13]) are available on Axiom A basic sets as they were for Anosov dynamical systems. 

4. SRB states on Axiom A attractors (see [92], [79], [14]). 

Let us cut an attractor B into a finite number of small cells such that the unstable 
manifolds are roughly parallel within a cell. Each cell is partitionned into a continuous 
family of pieces of local unstable manifolds, and we obtain thus a partition (S a ) of B into 
small pieces of unstable manifolds. If p is an invariant probability measure on B (for the 
dynamical system (/*)), and if its conditional measures a a with respect to the partition 
(Sq,) are a.c. with respect to the Riemann volume da on the unstable manifolds, p is called 
an SRB measure. 

The study of SRB measures is transformed by use of symbolic dynamics into a problem 
of statistical mechanics: one can characterize SRB states as Gibbs states with respect to 
a suitable interaction. Such Gibbs states can in turn be characterized by a variational 
principle*. In the end one has a variational principle for SRB measures, which we shall now 
describe. It is convenient to consider a general basic set B (not necessarily an attractor) 
and to define generalized SRB measures. First we need the concept of the (time) entropy 
/i/(/i), where / = f 1 is the time 1 map for our dynamical system, and p an /-invariant 
probability measure on B. This entropy (or Kolmogorov- Sinai invariant) has the physical 
meaning of mean information production per unit time by the dynamical system (/*) in 
the state p (see for instance [6] for an exact definition). The time entropy h/(p) should 
not be confused with the space entropy S and the entropy production rate e/ which we 
have discussed in Chapter 1. We also need the concept of expanding Jacobian J u . Since 
T x f maps E x linearly to EJ X , and volume elements are defined in E x , EJ X by a Riemann 
metric, we can define a volume expansion rate J u (x) > 0. It can be shown that the 
function log J u : B — > R is Holder continuous. We say that the (/^-invariant probability 
measure p is a generalized SRB measure if it makes maximum the function 

p ~ h f (p)-p(logJ u ) (9) 

One can show that there is precisely one generalized SRB measure on each basic set B; it 
is ergodic and reduces to the unique SRB measure when B is an attractor. The value of 
the maximum of (9) is precisely when B is an attractor, it is < otherwise. 

If m is a measure absolutely continuous with respect to the Riemann volume element 
on M, and if its density m vanishes outside of the basin of attraction of an attractor B, 
then 

p = w.lim^ i / dr(rym (10) 
t Jo 

* We shall come back to this topic in Chapter 3 (Section 3.4) 
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defines the unique SRB measure on B. We also have 

1 f* 

p = w.hm^oc - / drdf Tx 
1 Jo 

when x is in the domain of attraction of B and outside of a set of measure for the 
Riemann volume. The conditions (1-3) and also (6) are thus satisfied, and the SRB state 
p is a natural nonequilibrium steady state. Note that if B is a mixing* attractor, (10) can 
be replaced by the stronger result 

p = w.lim^oo (/*)*m 



5. SRB states without uniform hyperbolicity (see [66], [67], [95]). 

Remarkably, the theory of SRB states on Axiom A attractors extends to much more 
general situations. Consider a smooth dynamical system (/*) on the compact manifold 
M, without any hyperbolicity condition, and let p be an ergodic measure for this system. 
The theorem of Oseledec [69], [82] permits the definition of Lyapunov exponents X± < 
■ ■ ■ < AdimM which are the rates of expansion, p-almost everywhere, of vectors in TM. 
The Xi are real numbers, positive (expansion), negative (contraction), or zero (neutral 
case). Pesin theory [71], [72] allows the definition of stable and unstable manifolds p- 
almost everywhere. These are smooth manifolds; the dimension of the stable manifolds is 
the number of negative Lyapunov exponents while the dimension of the unstable manifold 
is the number of positive Lyapunov exponents. Consider now a family (S Q ) constituted of 
pieces of (local) unstable manifolds, and forming, up to a set of p-measure 0, a partition 
of M. As in Section 4 above we define the conditional measures a a of p with respect to 
(T> a ). If the measures a a are absolutely continuous with respect to the Riemann volumes 
of the corresponding unstable manifolds, we say that p is an SRB measure. 

For a C 2 diffeomorphism, the above definition of SRB measures is equivalent to the 
following condition (known as Pesin formula) for an ergodic measure p: 

h(p) = sum of the positive Lyapunov exponents for p 

(see [66], [67] for the nontrivial proof). This is an extension of the result of Section 4, 
where the sum of the positive Lyapunov exponents is equal to p(log J u ). Note that in 
general h(p) is < the sum of the positive exponents for the ergodic measure p (see [80]). 

Suppose that the time t is discrete (diffeomorphism case), and that the Lyapunov 
exponents of the SRB state p are all different from zero: this is a weak (nonuniform) 
hyperbolicity condition. In this situation, there is a measurable set S C M with positive 
Riemann volume such that 

^ n—l 

lim -V(5 fl = p 

n— >oo 77, z — ' 

k=0 

* i.e., for any two nonempty open sets Oi, O2 C B, there is T such that 0\ fl f l 02 ^ 
for \t\ > T. 
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for all x e S. This result [75] shows that p has the properties required of a nonequilibrium 
steady state. One expects that for continuous time (flow case), if supp p is not reduced to 
a point and if the Lyapunov exponents of the SRB state p except one* are different from 
0, there is a measurable set S C M with positive Riemann volume such that 



when x G S. I do not know a proof for this, however. 

After Pesin's theory, a more geometric approach to dynamical systems without uni- 
form hyperbolicity has started to develop (work of Jakobson, Benedicks, Carleson, Young). 
This is a difficult topic, which will not be dealt with here, but we refer the reader to the 
excellent review of L.-S. Young [99]. The interest of the geometric approach is that it gives 
more detailed results than the general Pesin theory: the existence of an SRB measure is 
proved instead of being assumed and the (usually exponential) decay of correlations can 
be studied, see Young [97], [98]. There is currently a rapidly growing literature on this 
topic, see for instance Viana [95], Chernov [15]. 

6. Rigorous versus formal arguments. 

We have just outlined some possible mathematical descriptions of nonequilibrium 
steady states. In the present situation however, where the theory is still in its infancy, 
priority should be given to getting the physics right. For the mathematical treatment, two 
main attitudes appear reasonable: 

(a) We may make uniform hyperbolicity assumptions - that the time evolution is 
given by an Anosov or Axiom A diffeomorphism or flow - and proceed rigorously. Of 
course uniform hyperbolicity is quite unrealistic from a physical point of view. 

(b) We may proceed formally and assume when needed that limits exist, that Lya- 
punov exponents are nonzero, or that correlation functions decrease fast. 

There are several reasons not to press too hard now the mathematical study of nonequi- 
librium steady states in the general situation of Section 5. One is that the mathematics of 
nonuniformly hyperbolic systems is still under development (for an important recent result 
see [5]). Another is that we shall, for physical reasons, want to take derivatives of SRB 
states with respect to parameters (Chapter 4), while these derivatives appear to make no 
mathematical sense in general. The resolution of this difficulty may involve the need to 
take a large system limit (thermodynamic limit) but this remains to be investigated. 



There is one zero exponent corresponding to the flow direction. 




X 
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Chapter 3. THE GALLAVOTTI-COHEN FLUCTUATION THOREM. 



1. Generalities. 

In this Chapter we prove an asymptotic formula due to Gallavotti and Cohen [44], 
[45] for fluctuations of the entropy production . The success of this formula supports 
the philosophy that calculations based on very strong hyperbolicity assumptions are in 
agreement with experiments (in this case numerical experiments). Following Gallavotti 
and Cohen we take the time to be discrete (t G Z), and the evolution to be given by an 
Anosov diffeomorphism (see Section 2.2) on the compact connected manifold M. We shall 
follow the ideas of Gallavotti and Cohen but change slightly the presentation. Instead of 
Markov partitions and symbolic dynamics we shall use a more geometric approach based 
on shadowing (Bowen [13]). 

2. Expansiveness, a-pseudoorbits, shadowing, and specification for Anosov 
diff eomor phisms . 

An Anosov diffeomorphism is expansive, i.e., there is 5 > (called expansive constant) 
such that 

((VfceZ) d(f k xj h y)<5) ^ x = y 

[In view of (2.8), when two orbits are close to each other, they must separate exponentially 
in the past or in the future, unless they coincide]. In fact we have, for some K > 0, 9 > 1, 

(d(f k x,f k y) < 6 when < n) d(x,y) < K8~ n (1) 

Let S be an interval of Z, finite or not, and let x = (xk)kes £ M s '. Given a > 0, we 
say that x is an a-pseudoorbit if 

d(fxk, %k+i) < ot whenever k, k + 1 G S 

We say that the orbit of x G M a- shadows x if 

d(f h x 7 x k ) < a for all k G S 

Lemma. (Bowen's shadowing lemma [13]). 

Given (3 > there is a > such that every a-pseudoorbit is ^-shadowed by the orbit 
of a point x G M. 

Note that if the pseudoorbit is periodic, then the true orbit which shadows it over Z 
is unique by (1), and periodic with the same period. [We assume that 2(3 is an expansive 
constant] . 

Given a > there is an integer R > such that for every pair (x, y) of points of M, 
and n > R, there is z such that 

d(z, x) < a , d(f n z, y) < a 
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This follows from the fact that / is topologically mixing (because we assumed / Anosov 
and M connected, see Section 2.2), and says that there is an orbit which interpolates in 
time n between any two points x and y. Using the shadowing lemma, we see that there is 
a true orbit interpolating between any number of a-pseudoorbits specified on on intervals 
with distances > R. This interpolation property is known as specification (see Bowen [10]). 

Lemma. (Specification property). 

Let a-pseudoorbits Xj be given on disjoint intervals of Z with mutual distances at 
least R. Then the x^ are all f3 -shadowed by the orbit of the same point x E M . 

3. Approximation of sums along orbits. 

We shall later have to consider sums of the form 

b 

$>(/ fc *) 

k=a 

If f k y stays close to f k x for a < k < 6, we know by (1) that d(f k x, f k y) < K9~ n where 
n = min(/c — a, b — k). Therefore, if A is Holder continuous* 

b b 

\J2Mf k x)-J2 A (f k y)\< L 

k=a k=a 

where L depends only on A. 

Suppose now that we use specification to obtain from x and y a periodic point z with 

f n z = z and 

f k z close to f k x for k = 0, . . . , r - 1 
f k z close to f k y for k = r + R - 1, . . . , n - R 
then we also have 

n— 1 t—1 n—R 

|$>(A)-]Ta(Ao- Mf k v)\<L' (2) 

fc = fc = fc = T + i?-l 

where V depends only on A. 

4. Periodic orbit representation of Gibbs states and SRB states for Anosov 
diffeomorphisms . 

We assume from now on that / is (at least) of class C 1+a for some a > 0. If A : M — > R 
is Holder continuous it is known that there is a unique /-ergodic measure p making 



h(p) + J A(x)p{dx) 



* Remember that A is 7-Holder continuous (with < 7 < 1) if \A(x) — A(y)\ < 
Cd(x,yr. Thus \A{f k x) - A{f k y)\ < Cd(f k xJ k V r < Ci^)-. 
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maximum; p is called the equilibrium state or Gibbs state for A (see Ruelle [81]). 

We need at this point the unstable Jacobian J u (see Section 2.4). By definition 
(rp^/\diraE u ma p S the volume element of E% to J u (x) times the volume element of Ej x , 
where the volume elements are defined with respect to the Riemann metric, and J u (x) > 0. 
Our assumption that / is C 1+a implies that J u is Holder continuous. In the present 
situation of an Anosov diffeomorphism, the SRB state was shown by Sinai [92] to be the 
Gibbs state pf for A = — log J u . Note that we have (for historical reasons) a somewhat 
confusing terminology: nonequilibrium steady states for a certain problem are described by 
equilibrium states (or Gibbs states) for another (very different) problem. There is however 
no ambiguity when we speak of SRB states. 

Let us write Fix/ n = {x G M : f n x = x}. An f-invariant probability measure cr n , 
carried by periodic orbits, is defined by 

^ ExeFix/- ®( x ) ex P Y2=l A (f k %) 
a n{^) = — rrrr^ 



Ex 6 Fix/™ eX P 

and it is known that, when n — > oo, a n tends weakly to the Gibbs state for A {i.e., cr n ($) 
tends to p(<E>) for all continuous $, see [81]). In particular, the SRB state pf satisfies 



ExeFix/^m^o 1 ^ 



X 



p f m = lim ±^±™ v /x " =u — (3) 

Z^xEFix/™ llfc=0 J \J X ) 

5. Reversibility. 

We say that the time evolution (/*), with t G Z, or R is reversible if there is a 
diffeomorphism i : M —> M such that i o i = identity and to/' = f~ l o i. [Note that 
any i o f s then satisfies the same conditions: tofojof = iojo f~ s o f s = identity, 
and i o / s o / = j o / s+1 = / _1 ojof]. We choose, as we may, a Riemann metric on M 
invariant under i. This gives in particular 

J(ifx) = J(x)- 1 (4) 

and, since i interchanges the stable and unstable directions, also 

J s (ifx) = J u (x)- 1 (5) 

For a system of particles with positions Xk and momenta pk one may take i((xk), (Pk)) = 
((xk),(—pk)) if the forces are not velocity dependent (electrostatic forces are permitted, 
but not magnetic fields). In this physical situation one speaks of microscopic reversibility. 
This was an essential ingredient in Onsager's proof of his reciprocity relations, and it is 
also essential for the fluctuation theorem of Gallavotti and Cohen. 

6. Entropy production. 
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We have discussed earlier (Section 1.9) the formula 



e/ = -/ W (i)lo g J(x)>0 



for the average entropy production per unit time in the nonequilibrium steady state defined 
by the SRB measure pf. Remember that J > is the Jacobian of / with respect to the 
chosen metric on M. In the situation which interests us here, where / is an Anosov 
diffeomorphism of the connected manifold M, we have e/ > unless / has an invariant 
measure absolutely continuous with respect to the Riemann volume. [This follows from 
Ledrappier [65], see the discussion in [85]]. 

For integer r > 0, let 

e T {x) = —Y j \ogJ{f k x)- 1 (6) 

7 k=0 

This may be interpreted as a dimensionless entropy production rate or phase space con- 
traction rate at x over time r, such that p/(e T ) = 1. Note that if we have reversibility, (4) 
yields 

e r o % o f = -e T (7) 

The fluctuation theorem of Gallavotti and Cohen (Theorem 9 below) is a statement about 
the fluctuations of e r with respect to the SRB measure p / (or equivalently with respect to 
time). 

7. Proposition. 

Let f be a C 1+Ct Anosov diffeomorphism of the compact connected manifold M , and let 
Pf be the corresponding SRB measure. Assume reversibility (io f = f~ x oi, i 2 = identity ) 
and define e T by (6) with respect to an i-invariant Riemann metric on M. Under these 
conditions there are constants a, b > such that for any interval [p, q] C R, 

1 log P/({x:e T (x)e[p, Q ]}) ^ q+ b (g) 



re/ p f ({x : e T (x) G [-q-a/r, -p+a/r}}) t 

Let x De the characteristic function of [p, q] and x* the characteristic function of 
[p — a/r, q + a/r] where a will be fixed later. Using (7) we have 

p f ({x : e T (x) G \p,q}} _ p f ( X oe T ) 



p f ({x : e T (x) G [-q - a/r, -p + a/r]} p f (x* o (-e T )) 

= P/(X°gr) = P+(X °er) 

P/(x*°er°«) p_(x*oe T ) 

where p + = pf and p- = i*pf is the anti-SRB state, i.e., the Gibbs state for log J s (where 
J s is the stable Jacobian). We replace now x, x* by continuous functions 0, 0* (we shall 
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take (p I x, (p* I x* later) and use the formula (3) (approximation of SRB states by periodic 
orbits) to rewrite 

P+(0°Ct) 

P-{4>* °er) 

ExeFix/" 0( e r(^)) Ilfc=o J u (f k x)~ 1 E x eF^<P*(cAx))ll n k ZoJ s (f k x) 



lim 



E^eFix/™ rifc=o J u (f kx ) 1 ExeFix/" rifc=o J s (f kx ) 

= lim 



E,6Fix/" 0*(^(^))nfe=O^(/^) 

[The last step is based on (5) and the fact that i is a bijection Fix/ n — * Fix/ n which maps 
orbits of / to orbits of / (reversing the direction)]. 

If x G Fix/ n , let S(x) C Fix/ n consist of those y such that f k y is close to f k x 
for k = 0, . . . , r — 1, and to if n ~ k+T x for k = r + R — 1, . . . , n — R. By specification 
cax&S(x) > 1, and by expansiveness cax&S(x) < N for some N independent of n, r. 
Suppose y G S(x), in view of (2) and (5) we have 

n— 1 t — 1 n—1 

n ^(/ fc 2/) _i < m n n 

k=0 k=0 k=r 

with M independent of n, r. We also have 

a 



\zAv) ~ *A X )\ < 



T 



with a independent of n, r (this fixes our choice of a). Therefore x{ e Av)) — X*( e r( x )), 
and we may assume that the continuous functions 0, </>* (such that 4> I Xi 4>* i X*) a l so 
satisfy 0(e T (y)) < 0*(e T (x)). Thus 

n-1 

xeFix/™ fe=0 

n-1 

< E E ^Ay^U^ify)- 1 

xeFix/™ y£S(x) fc=0 

r— 1 n—1 

< E <^(x)<f>*(e r (x))M]lJ»(f k x)- 1 ]lJ s (f k x) 

xEFix/™ fc=0 fc=r 

r— 1 n—1 

<mn e 0*(e.^))n ju (/^) _i n js (/ fe:r ) 

xeFix/™ fc=0 fe=r 
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Hence 

-™ E« e EW-^(Cr(x))lUi;j-(/ fc x)lK^J-(/*x) 

We still need the fact that 



> p 



where the constant P > is independent of x and r (this follows from the continuity of 
the splitting TM = E s ©E", and compactness). Thus, by (6), 



< -ex.p(re f e T (x)) 



Inserting in (9) and performing the limit n^oowe find, using Section 4 and the definition 

of p_, 

p+((j)oe T ) MArp_((0*oe r ).exp(re/e T )) 



p_(</>*oe T ) P P-{4>* ° e r) 

Letting now 0* j x* and using e T (x) < q + a/r when x*(e T (x)) 7^ 0, we obtain 

P+(x°tr) ,MN a 

— r~* V - ~5~ exp re/ g + -)) 

p_(X*oe T ) P r 

which yields (8) if we take b = a + (l/e f ) \og(MN/P). □ 
8. A large deviation result. 

Proposition 7 is a first form of the fluctuation theorem. To obtain the final form we 
need the following large deviation result. 

Let A, B : M — > R be Holder continuous, and let p be a Gibbs state for A. There 
exists then a real concave function 77 on the open interval (pi, £><■>) sucn that, for every 
interval I with I n {p\,pl) 7^ 0, 

1 1 r_1 

lim -log p({x : B(f h x) e I}) = sup 77(14) (10) 

T ^°° T T t^o ueJn(p; )P 3) 

(We disregard here a degenerate case where pi = p\ and (1/r) X]fc=o B(f k x) does not 
fluctuate for large r]. 

By use of symbolic dynamics, this is converted into a fact of equilibrium statistical 
mechanics, for which see Lanford [59]. Because we have here one-dimensional statistical 
mechanics with short range interactions, 77 is in fact real analytic and strictly concave. 
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We take now A = — log J u so that p = pf, and we let B = — (1/e/) log J. We have 

thus (1/t) YlkZo B(f k x) = e T {x). Also, we have (p*,^) = ( — lAP*)? *.e., e r fluctuates in 
an interval symmetric around 0; this follows from reversibility as should become clear in a 
moment. Equation (10) now has the form 

lim \ogpf({x : e T (x) G /}) = sup fj{u) 

r^oo re/ uein(- P *, P *) 

(with fj = rj/ef). In particular, if ]p| < p* and 5 > we have 

lim — log/9f({x : e T (x) G (— p — 5, — p + 5)}) 
r^oo re/ 

= lim — logp/({x : e T (x) G (-p - 5 - -, -p + S + -)}) 
t^oo re/ r r 

From (8) we obtain also 

lim J_ log Pf({x:e T (x)e(p-S,p + S)}) ^ 

r^oo re/ P/({z : e T (x) G (— p — 5 — a/r, — p + 5 + a/r)}) 

Therefore 

lim J_ log />/({* = M*)e fez^£±j|L < p + i (11) 

r^cx) re/ P/({^ : e T (x) G (— p - o, -p + b)\) 

9. Theorem (Gallavotti- Cohen fluctuation theorem [44], [45], [36]). 

Let f be a C 1+a Anosov diffeomorphism of the compact connected manifold M , and let 
Pf be the corresponding SRB measure. Assume reversibility (io f = f~ x oi, i 2 = identity ) 
and define e T by (6) with respect to an i-invariant Riemann metric on M. Then there 
exists p* > such that if \p\ < p* and 5 > 

c- ^ t It Pf (i x '• e r{x) G (p - 5,p + 5)}) 
r^cvref Pf{{x : e T {x) G (— p - o, -p + o)\) 



This follows from (11) and the formula 

lim J_ lQg Pfii^ ; gr(x) G (-p - - p + fl}) + ^ 

r^ooref p f ({x : e T (x) E (p - 5,p + 5)}) 

obtained by changing p to — p in (11). [j 

10. Extension of the fluctuation theorem to time-reversed events. 

Inspection of the proof of the Gallavotti-Cohen fluctuation theorem shows that, basi- 
cally, it compares the probability P + of the event e T {x) xs p with the probability P_ of the 
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time reversed event e T (if T x) pa p, obtaining P + /P_ pa e e ^ pr . As remarked by Gallavotti 
[43], one might as well compare (for any function W) the probability P + of the event 

e T (x) pa p and \I/(:r) « ^4 

with the probability P_ of the event 

e T (if T x) pa p and ^(if T x) pa A 

again with the result P + /P- ~ e efpT . Note that ~ A could for instance describe 

density fluctuations of a system of particles in the time interval [0,r], and ^ oio f T (x) pa A 
the time reversed fluctuations in the same interval [0,r]. Also note that the right hand 
side of P + /P- pa e efpT depends on p but not on A. A specific example of an extension 
of the Gallavotti- Cohen fluctuation theorem is the following formula, where \1> is a strictly 
positive continuous function on M: 

1 ft iw Pf(dx)ty(x) 

6 < lim \ J l e -r-p|<o rj v _____ < _|_ ^ 



T- 



Remember that p/(e r ) = 1, so that e T takes mostly positive values. Gallavotti [43] 
points out that in those exceptional cases when e T (x) ~ — p < one observes in the 
interval [0, r— 1] (roughly speaking) the time reversed behavior from the "normal" behavior 
corresponding to e T (x) ~ p < 0. This means that the pattern of time fluctuations is 
reversed, i.e., the time appears to flow backwards. 

11. Other extensions. 

Note that the results of this chapter apply far from equilibrium (we assume that 
the time evolution is reversible, but we do not assume that it is close to Hamiltonian) . 
Placing ourselves close to equilibrium we could obtain from the fluctuation theorem a 
linear response formula called fluctuation- dissipation theorem, and the Onsager reciprocity 
relations (see Gallavotti [37]). We prefer to postpone the derivation of these results, and 
obtain in the next chapter a linear response formula far from equilibrium and without 
assumption of reversibility. 

Naturally, one wants to extend the fluctuation theorem from Anosov diffeomorphisms 
to Anosov flows. Such an extension has been obtained by Gentile [49], using Markov parti- 
tions. Probably one could use instead an approximation of SRB states by periodic orbits, 
and Bowen's techniques (see [11] and [70]). 

Physically, the assumption of an Anosov time evolution is not as unrealistic as it would 
at first sight appear. In fact, close to equilibrium, hard sphere systems are expected to be 
hyperbolic (and differ from Anosov systems only by the presence of singularities). Problems 
arise, however, as we move away from equilibrium. In particular, time reversibility implies 
that the stable and unstable dimensions {djmE s and dimE u ) are equal to dimM/2 for an 
Anosov diffeomorphism (or (dimM — l)/2 for an Anosov flow). Therefore the numbers 
of positive and negative Lyapunov exponents are also equal (for any ergodic measure, in 
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particular the SRB measure pf). If pairing is satisfied (see Dettmann-Morriss [22]) the 
Lyapunov exponents come by pairs A_j, \ with A_i+A^ = 2a independent of i (for the SRB 
measure pf, we have a = — e//dimM, resp. a = — e//(dimM — 1) for a diffeomorphism 
resp flow). The Anosov assumption thus forces a gap of size 2\a\ between the positive 
and the negative exponents, which seems unphysical. To avoid this problem, Bonetto and 
Gallavotti [7] consider a situation where / satisfies the Axiom A of Smale [93] , but the SRB 
measure pf and the anti-SRB measure i*pf may now have disjoint supports. Assuming a 
new symmetry (obtained from what is called Axiom C in [7]), Bonetto and Gallavotti can 
still derive the fluctuation theorem. In the absence of uniform hyperbolicity, the number 
of positive Lyapunov exponents may vary discontinuously, and it is difficult to relate pf 
and i* pf even if they have the same support; a fluctuation theorem has not been derived 
under these more general circumstances. 

The available numerical simulations (by Evans, Cohen, Morriss [31], and Bonetto, 
Gallavotti and Garrido [8]) are in excellent agreement with the fluctuation theorem, but 
in a situation where the equality between the numbers of positive and negative Lyapunov 
exponents is respected. Note that, in the form given above, the fluctuation theorem ap- 
plies to small systems; for large systems, the entropy production fluctuates very little, 
and negative values are unobservably rare. Following a suggestion by Gallavotti, there is 
currently a lot of interest in a possible local fluctuation theorem (which would apply to 
small subsystems of a large system). For other interesting recent work, see Kurchan [58], 
Lebowitz and Spohn [64]. 
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Chapter 4. LINEAR RESPONSE. 



1. Setting the problem. 

We shall see that the definition of SRB state extends to time dependent dynamics 
(nonautonomous systems). This gives in particular the possibility of studying the re- 
sponse of a stationary nonequilibrium state to small time dependent perturbations. The 
linear response function thus defined is a physically important quantity, and its study is 
the main subject of this chapter. Our approach to computing the linear response will be 
mostly formal. Actually, in the uniformly hyperbolic case, one can rigorously establish the 
differentiability implied by linear response, but it should be realized that this situation is 
exceptional. In general, i.e., unless we are in some kind of structurally stable situation, 
changing the parameters of a dynamical system (unless the changes are in special direc- 
tions) will entail rapid discontinuous changes in the dynamics, and also in the SRB states. 
From a physical point of view this may not be as bad as it seems, because in physical 
applications we are mostly interested in large systems. For such systems (i.e., in the ther- 
modynamic limit) one can hope that the dependence on parameters is in general smooth, 
and that our linear response formula is physically correct. To put things more bluntly, if 
there is a linear response formula at the level of experimentally observable quantities, it 
must be (10) as discussed below. This is because our derivation, altough formal, offers no 
reasonable alternative. What may happen is that the integral in the right-hand side of 
(10) diverges, leaving us without a response formula. 

The discussion of the present chapter is valid far from equilibrium, i.e., valid in general. 
The assumption of closeness to equilibrium and its consequences will be considered in the 
next chapter. While our derivation of the linear response formula will be formal, and easy, 
it involves a nontrivial element, viz., the SRB state p. By contrast with the situation close 
to equilibrium, where one can make guesses on the basis of general physical principles, we 
shall in what follows need specific properties of SRB measures. 

2. Time dependent SRB states. 

It will be convenient to assume for the purposes of this section that we have a discrete 
time t e Z. Instead of a single diffeomorphism / of the compact manifold M we consider 
a family (ft)tez, i-e., we have time dependent dynamics. What we want to call SRB state 
in this situation is a family (pt)tez where p t is a probability measure on M, such that 

p t = lim f* ... ft_ n m (1) 

n^oo 

and the probability measure m is absolutely continuous with respect to the volume element 
on M. In particular p t = flp t _\. 

Note that the above situation arises for random dynamical systems* where a shift- 
invariant measure P(da) is given on the a = (ft)tez, and statements are made for P-almost 

* Such random dynamical systems have been much studied, see Kifer [55], Arnold [2], 
and for the SRB states Bahnmiiller and Liu [3]. 
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every a. (The case of a perturbation proportional to e wt that we shall consider later can 
be handled in this setup). 

For applications to physics, however, it is desirable to be able to consider a nonprob- 
abilistic situation. Secifically, we proceed with a brief discussion of the case where all ft, 
t e Z, are sufficiently close to an Anosov diffeomorphism* / of M. Structural stability 
extends to this situation and gives a (unique) family (jt)tez of homeomorphisms of M 
close to the identity such that we have a commutative diagram 

••• -> M h M M -> ••• 

T jt-i T n T h+i 

••• -> M X M M — > • • • 

Using hyperbolicity, the fact that the jt are close to the identity, and Tf t close to Tf, one 
can define stable and unstable subbundles V t s , of TM with the usual properties. The 
corresponding (local) stable and unstable manifolds V t s (x), V"(x) are smooth because we 
are in a uniformly hyperbolic situation. One checks also readily that jt~ 1 Vt' u (jtx) coincides 
(in a sufficiently small neighborhood of x) with the stable or unstable manifold V s,u (x) for 
/. Note that although j t is not smooth in general, it maps the smooth manifolds V s,u {x) 
to the smooth manifolds Vt' u (x). Let us now define the time dependent SRB state (p t ) 
by (1). One can then show that the conditional measures of pt on the unstable manifolds 
V™(x) are absolutely continuous with respect to the Riemann volume element of V"(a?) and 
satisfy a uniform continuity condition in t. Conversely, if (ptjtez satisfies flp t -\ = pt and 
the above absolute continuity and uniformity, then (p t ) is a time dependent SRB state. 
Such states can also be studied by transfer operator techniques (see [86]). 

Although time dependent SRB states are conceptually important for nonequilibrium 
statistical mechanics, we shall mostly use those which correspond to an infinitesimal time 
dependent perturbation of a time independent dynamical system. Also, we shall mostly 
use a continuous rather than discrete time, and work formally rather than rigorously. 

3. A general linear response formula. 

Consider time dependent dynamics defined by 

% = F ^ < 2 > 

on the compact manifold M. Integrating (2) with initial condition x at time s gives at 
time t a point x(t, s) = f(t, s)x. For this time evolution, a time dependent SRB state 
(pt)teR is defined by 

p t = lim f(t, s)*m (3) 

s—> — oo 

where m is absolutely continuous with respect to the volume element on M, and we assume 
that the limit (3) exists. 



* More generally one could consider an Axiom A diffeomorphism / near an attractor, 
see [86] , Section 4 
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Let now (p t + 5p t )te~R. °e the time dependent SRB state corresponding to the time 
evolution 

- = F t (x) + 8F t (x) 

we have thus 

p t (<f>) + 8p t (<f>) = lim m(Qo(f(t,s) + 8f(t,s))) 

s^ — oo 

Proceeding formally this gives 

5pt(&) = lim / m(dx)(V x(ttS) $) ■ Sx(t, s) 

s— > — oo J v ' 

and we may assume s < t. Using the tangent map T x f to / at x we may write 

Sx(t,s)= / dT(T x ( Tt8 )f(t,T))6F T (x(T,s)) 

J s 



hence 



6pt($) = dr j m(dx)((T x{TiS) f(t,T))6F T (x(T,s)))-V x{t , s) $ 

= s hm oo j'dr J(f(r, s)* m) (dy) ((T y f (t , r))SF T (y)) ■ V B(t , T )$ 
dr J p T (dy)((T y f(t,T))6F T (y))-V y{t , T) $ 



t 



Finally 

6 t p(*) = [ dr j p T (dy)SF T (y) ■ V y ($ o /(*, r)) (4) 

J — oo J 

= [ dr j p t (dx)((T xiTjt) f(t,T))6F T (x(T,t))-V x $ (5) 

J —oo J 

This is the most general form of our linear response formula. 

The convergence of the right-hand side of (4), (5) is assured if r — > 8F T has compact 
support, but one would like to assume just boundedness, to allow time independent and 
time periodic perturbations. As discussed above, one expects that (4), (5) can be proved 
under uniform hyperbolicity assumptions, but not in general. As things stand, the case of 
Anosov diffeomorphisms and time independent perturbations has been treated rigorously, 
and one obtains the following result. 

4. Theorem (Ruelle [86]). 

Let A be the space of C r Anosov diffeomorphisms of the compact connected manifold 
M, and let p be the SRB state for f e A. If r > 3, the map f i-> p\C r ~ 1 (M) is 
C r ~ 2 : A — > C r_1 (M)*. 
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The derivative of / i— > p is given by 

oo 

5#) = ^KX'V($of)) (6) 



n=0 



where we have introduced the vector Geld X = 5f o / 1 , and taken $GC r 1 (M). 

Note that for time dependent dynamics (ft)tez and time dependent SRB state (p t ), 
(6) should be replaced by* 

t 

M*) = E Pr(*T-V(*0/ t 0...0/ T+1 )) (7) 

r= — oo 

which is the discrete version of (4). 

The existing proof of Theorem 4 applies actually to Axiom A attractors, and even to 
general "Axiom A basic sets" , provided the definition of SRB measures is generalized and 
(6) replaced by a more complicated formula (see [86]). 

5. Time dependent perturbation of a nonequilibrium steady state. 

We shall now assume that F t = F is time independent, and let pt = p be a time 
independent measure corresponding to the time evolution 

We may here write /(£, s) = / t_s , x(t, s) = x(t — s) = f f ~ s x. The pertubation 5F t , noted 
now S t F, will remain time dependent, so that the perturbed time evolution is 

^ = F(x) + 6 t F(x) (9) 

Let (p + 8 t p) be the time dependent SRB state for the evolution (9) which replaces the 
time independent SRB state p for the evolution (8). We are interested in the response 
(5 t p) to the perturbation (StF). In view of (4), (5) we have (see [88]) 

M*) = / dr f p(dy)S T F(y) ■ V y ($ o f~ T ) (10) 

J — oo J 

= f dr j p(dx)((T x{T _ t) f- T )5 T F(x(T-t))-V x $> (11) 

J — oo J 

Assuming that $ is in a suitable space B of functions on M, and SF in a suitable space X 
of vector fields we may rewrite (10), (11) as 



StP 



= J dr K t - T 5 T F = J da K a 5 



* At the time of writing, there is no written proof of (7). 
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where the linear operator n a maps X to the dual B* of £>, and 

(K a X)$ = for a < 

(K a X)$ = J p(dx)((T x{ _ a) r)X(f-°x)) ■ V x $ for a > (12) 

6. Properties of the response function. 

The (operator-valued) function a i— > k ct defined by (12) is called response function. 
The fact that K ff = for a < is called causality and is interpreted to mean that the cause 
8 T F must precede the effect 8 t p. 

To study the behavior of K a for large a we shall assume a weak form of hyperbolicity, 
namely that all Lyapunov exponents* with respect to p are ^ 0, except one corresponding 
to the direction F of the flow. We may thus write for any vector field X on M 

X = X s + (f>F + X u (13) 

where X s is in the stable (contracting), X u in the unstable (expanding) direction and <j> 
is a scalar function. 

Lemma. (Ruelle [88]). 

For a > we have 

(a)$ = J p(dy)[((T y nX s (y))-V y{a ^-(F-V<P + dw u X u )(yMy(a))} (14) 

where the unstable divergence div u will be defined below. 
We first obtain from (12), (13), 

M)* = / p(dy)[((T y nX s (y)).V y{a ^ + ^F + X u )(y).V y ^or)] 

To transform the </>-term we use 

<j>F ■ V* = <t>jV o /*| t=0 = ^((0*) o f% =0 - (^0 o /*| t=0 )tt 

= ^((0*)o/*)| t=o -(F.V0)« 
so that with \1> = $ o / CT we get 

| p{dy){<j>F){y)-V y {®or) = - J p(dy)(F-Vct>)(y)3>(y(o-)) 



* See Section 2.5. 
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To transform the X"-term we shall use the fact that p is an SRB measure. We may 
write for p a disintegration of the form 



p{dy) = j u(d\)p {X) {dy) 

where p/ X )(dy) is a conditional measure on a local unstable manifold. The SRB condition is 
that p^x) (dy) is absolutely continuous with respect to the Riemann volume element d? x ^ (dy) 
on the (local) unstable manifold labelled by A. Writing p^(dy) = m^\)(y)d^(dy) one 
sees that the density m^(-) is unique up to a multiplicative constant, in fact 

m(x) _ A J u (f~ n x) 

~ J- = \ J u (f~ n y) 

where J u denotes the unstable Jacobian* for /. Let us place ourselves in the unstable 
manifold indexed by A. Since p( X ) (dy) is a volume element, its Lie derivative with respect 
to a vector field Y satisfies 

LyP(\) = (dw u Y)p {x) 

where div u y is the divergence of Y with respect to the volume element p^(dy). Note that 
div u does not depend on the choice of arbitrary multiplicative constant in the definition 
of p(\^(dy). We call div u the unstable divergence and obtain 

J p W (dy)(Y-V*)(y) = - J p (x) (dy)(dw u Y)(y)^(y) + boundary term 
Since X u is tangent to the unstable manifolds we have 

j p(dy)X u (y)-V y V = J i/(dA) J p {x) (dy)X u (y) ■ V w ¥ 

= - J v(d\) j p (x) (dy)(dw u X u )(y)^(y) = - j p(dy)(div u X u )(y)^(y) 

where the boundary terms disappear because of cancellation between the local pieces of a 
global unstable manifold. Taking \1/ = $ o f a we obtain thus 

p(dy)X u (y) • V y ($ o f-) = - J p(dy)(div u X u )(y)<S>(y(<j)) 

which completes the proof of (14). [! 

7. Physical discussion: the question of integrability of the response func- 
tion. 



* See Section 2.4. 
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We interrupt here for a moment our technical discussion, and we pause to consider 
the relation of our analysis of linear response with more traditional approaches. (For the 
latter see for instance de Groot and Mazur [20]). 

The usual approaches usually suppose closeness to equilibrium (as discussed in Chap- 
ter 5 below) and assume that the response of a system to small perturbations is empirically 
given. Physical principles such as causality or microscopic reversibility are then seen to 
have consequences for the response, and these consequences are investigated. 

In the analysis presented here, by contrast, the response of a system to perturbations 
is deduced from the general principle that nonequilibrium steady states are given by SRB 
measures. This new approach, being much more explicit, applies directly far from equilib- 
rium, and causality is a consequence of the choice of SRB measures instead of appearing 
as an extra assumption. 

In the traditional approaches it appears reasonable to assume that the empirically 
defined response function is well behaved, say integrable. In the new approach this has 
to be discussed on the basis of expressions such as (12) or (14), and we must admit that 
we are on somewhat shaky ground. Formally, we can argue as follows. Since X s is in the 
contracting direction, the term 



is a correlation function (note that J p(dy)(F-\7(j)+dw u X u )(y) = 0), and we may hope that 
it decreases exponentially when a — > oo. Therefore (K a X)3> is an exponentially decreasing 
function of cr, and is in particular integrable. 

In fact difficulties arise because the terms X s , <fiF, X u of the decomposition (13) of 
X need not be continuous functions on M unless we assume uniform hyperbolicity, and 
because even for Anosov flows the exponential decrease of correlations has been proved only 
under certain conditions (see Dolgopyat [25], [26]). That these difficulties are not visible 
in the traditional approaches is because one is satisfied there with an empirical definition 
of the response function. As soon as one tries to obtain a dynamical understanding of the 
response, difficulties necessarily arise. A resolution of these difficulties will require a better 
understanding of smooth dynamics, and a better understanding of the infinite system limit 
(thermodynamic limit). 

While we do not underestimate the difficulties just discussed, we shall want to push 
ahead, assuming that the response function may somehow be considered as integrable, as 
is done without trepidation in the more traditional approaches. The physical interest of 
the results obtained may then justify further concern for their mathematical justification. 

8. Dispersion relations. 




should decrease exponentially with a. The term 
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The Fourier transform of the response function is called the susceptibility: 

k w = / K a e iua da 



Remember that k u , k^ are linear operators X — > £>*. Choosing X G A 1 and $ G £>, we 
note that (k^X)^ is the response of a system in the nonequilibrium steady state p to a 
periodic perturbation. 

Taking a h- > (k ct X)$ to be integrable as discussed in Section 7, we see that cj h- > 
(£ ti ,X)<E> is continuous and tends to at infinity (Riemann-Lebesgue lemma). Furthermore, 
the fact that n a = for a < implies that a; i— > (£ ti) X)$ extends to an analytic function 
in the upper half complex plane. The analyticity of u> i— > {k u X)Q/{oj — it) and standard 
manipulations (contour integration) yield then 

(k u X)$ = — p.v. / —doj 

m J-oo uj-u 

where p.v. denotes principal value. More formally 

k u = — p.v. / — — du (15) 

TTl J_ 00 UJ-U 

Splitting k into real and imaginary parts, one can also write that one is the Hilbert trans- 
form of the other. The functional identity (15) is the Kramers- Kronig dispersion relation. 
Note that no condition of closeness to equilibrium was necessary to derive (15). 

9. Correlation functions, spectral density and resonances. 

Let p be any probability measure on the compact space M, invariant under the time 
evolution (f t )tez, or (/*)t e R. There is a group homomorphism U of Z or R into the 
unitary operators on L 2 (M, p), defined by 

U(t)V = #<)/* 

There is also a homomorphism a of the algebra C(M, C) of complex continuous functions 
on M into the algebra B(L 2 {M, p)) of bounded operators on L 2 (M, p), defined by 

(a(A)V)(x) = A(x)V(x) 

The homomorphism a satisfies a(A*) = (a A)* (adjoint operator) and a(l) = 1 (unit 
operator). Furthermore 

a{Aof) = U(t)a(A)U(-t) 
and if O is the function 1 considered as element of L 2 (M, p) we have* 

U(t)Q = n 



* The construction just indicated is a special case of the GNS construction (Gel'fand- 
Naimark-Segal). 
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(n,a(A)fl) =p{A) 
If A, B G C(M, C) we define the (time) correlation fuction 

p AB (t -s) = p((A o o /*)) - p{A)p{B) 

= (a(A*)Q, U(t - s)a(B)tt) - (fi, a(A)Q)(Q, a(B)Q) 
Introduce the spectral representation 




P(dw) 



where P is the spectral measure on the circle R(mod 2tt) in the discrete time case or on 
R in the continuous time case. Note that P({0}) is the orthogonal projection on the space 
Ho of U(t) -invariant vectors (for all r G N or R), and that O G Ho- Ergodicity of the 
measure p with respect to (/*) is equivalent to dim7io = 1 (i.e., there are no nontrivial 
invariant functions). Assuming ergodicity we have thus 

P({0}) = orthogonal projection Pq on f2 

Writing * x = a(A*)Q, # 2 = a(B)Q, we obtain 

Pab(t) = p(A(B o f T )) - p(A)p(B) = U(t)V 2 ) - P n * 2 ) 

J {u^O} 

In particular, if A = A* and \1> = a(A)Q we find 

Paa(t) = J e luJT paa(u) du 

where 

(*,P(dw)*) = (paa(u) + \(t>,n)\ 2 8(uj))daj 

and paa(u) > 0. The time autocorrelation function paa is thus the Fourier transform 
of the positive density paa(<^), called spectral density, or power spectrum (of the "signal" 
t^Aof). 

Let us now specialize to the situation where p is an SRB measure for a smooth dy- 
namical system (/*) on a compact manifold. Under various hyperbolicity assumptions 
(in particular uniform hyperbolicity) and under suitable smoothess assumptions on A, £?, 
one can show that pab is the Fourier transform of a continuous function 10 — > pab(w) 
which extends meromorphically to |Imo;| < c for some constant c > 0. The poles of pab 
are called resonances, their position is independent of the choice of A, B (although some 
residue might accidentally vanish) . For a discussion of resonances see in particular Ruelle 
[83], [84], Parry and Pollicott [70], Baladi [4], Dolgopyat and Pollicott [27]. 
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10. Fluctuation-Dissipation. 

In nonequilibrium statistical mechanics close to equilibrium there is an important 
result known as fluctuation-dissipation theorem (see Chapter 5). We discuss here briefly 
and formally what remains of this result far from equilibrium on the basis of the following 
lemma. 

Lemma. 

The susceptibility may be written as 
where 

/•OO 

= p[( / da{{TDX s ) o /-") ■ V*] 

Jo 

/•OO 

= - e lLja dap[((F ■ V</> + div u X u ) o 
Jo 

This follows directly from the lemma of Section 6. [] 

The functions to i— > (k*X)<E>, (k™X)$> extend analytically to the upper half-plane by 
causality. Formally, the function a i— > p[((F ■V4>+div u X u )o f~ a )^] is a correlation function 
(note that p(F ■ V0 + div^X") = 0); therefore we may hope that its Fourier transform 
extends meromorphically to a strip |Imo;| < c as discussed in Section 9. This implies that 
u> I— > (k™X)3> extends meromorphically to Imcu > — c, its poles for > Imo; > — c being 
the same resonances as those of the spectral density. 

Formally, we also obtain that (k^X)® is holomorphic for Imu > A where A is the 
largest negative characteristic exponent of (Tf a ) with respect to p. We can also define T a 
on the vector fields in the stable direction by 

(TrY)(x) = (T f -. x r)Y(f-'x) 

Then, on a suitable Banach space, (T cr ) cr > is a contraction semigroup and, if —H is its 
infinitesimal generator we have by th Hille-Yosida theorem 

/•OO 

/ e lLja da((Tf a )Y) o = (H - iuj)~ l Y 
Jo 

so that (k^X)® = p[((H — iu)~ 1 X s ) ■ V<&]. The singularities of are thus related to the 
spectrum of H. 

In conclusion the susceptibility has stable and unstable poles, which can in principle be 
distinguished experimentally because the unstable poles are the same as those of a spectral 
density. Close to equilibrium the difference between these two kinds of poles disappears 
as we shall see in the next chapter. 
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Chapter 5. NONEQUILIBRIUM NEAR EQUILIBRIUM. 



1. Terminology and assumptions. 

An isolated macroscopic system which is not subjected to external forces normally 
evolves towards a macroscopic state of rest* corresponding to certain values of thermody- 
namic variables like the temperature. At the macroscopic level, this state of thermody- 
namic equilibrium is described by the microcanonical ensemble of equilibrium statistical 
mechanics (see Section 1.9). Typically, we have a phase space R 2Ar with a Hamiltonian 
function (p, q) i— > H(p, q) where p = (pi, . . . ,Pn), q = (qi, ■ ■ ■ , Qn), and the so-called 
microcanonical ensemble corresponding to a certain energy K is the probability measure 

p(dpdc l ) = ^8(H(p,< l )-K) (1) 

where O is a normalizing factor and it is assumed (typically) that M = {(p, q) e R 2Ar : 
H(p, q) = K} is compact and that gradif does not vanish on M, so that M is a smooth 
compact manifold. 

Note that at the microscopic level we do not have rest or equilibrium, but a nontrivial 
time evolution 

1(P\- f-<Vf(p,q)\ 

dt{<i) \ <9 p #( P ,q) ) [Z) 

We shall nevertheless use the traditional terminology and call equilibrium the situation 
where we have a smooth time evolution (/*) (defined by (2)) on the compact manifold 
M, and a probability measure p (defined by (1)), invariant under (f*), and with smooth 
nonvanishing density with respect to the Riemann volume element dx on M. It is desir- 
able to assume that p is ergodic. In physical applications, p may have one large ergodic 
component, with the measure of the remainder tending to zero in the thermodynamic limit 
(i.e., for large systems); it would not be difficult to adapt our discussion to this situation. 
In the present chapter we shall disregard the symplectic structure implied by the Hamilto- 
nian time evolution (/*) (symplectic structure will be important for the Dettmann-Morriss 
pairing theorem, see Chapter 6). 

Notice that the equilibrium state p is an SRB measure. In fact, for the purposes of this 
chapter, an equilibrium state p may be defined as an SRB measure with smooth density 
0) with respect to dx, the time being allowed to be discrete (t e Z) or continuous (f e R). 
We shall see in Section 3 that the entropy production vanishes for an equilibrium state, 
which is therefore described as nondissipative. (An SRB state with entropy production > 
is dissipative) . In Section 7 we shall assume the reversibility of the dynamics (i.e., there is 
a diffeomorphism i : M — > M such that i 2 = identity, if* = f~*i), and of the equilibrium 
state p (i.e., i*p = p)**. Actually the latter condition follows from the assumed ergodicity: 

* We assume for simplicity that the total momentum and angular momentum vanish. 
** In Chapter 3 we assumed reversibility of the dynamics, but the "anti-SRB" state i* p 
was different from p. The condition i* p = p can be imposed only at equilibrium. 
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i*p is absolutely continuous with respect to dx, hence with respect to the ergodic measure 
p, and since i*p is (/*) -invariant, i*p = p. 

To summarize: We are given a smooth dynamical system (f f )tez or (f f )teR on the 
compact manifold M, with an ergodic probability measure p called equilibrium state. We 
assume that p(dx) has smooth nonvanishing density with respect to dx. In Section 7 we 
also assume reversibility, i.e., there is a diffeomorphism i of M such that i o i = identity, 
io/' = f~ l o i; in particular this implies i* p = p. 

These assumptions will allow us to recover the results of nonequilibrium statistical 
mechanics close to equilibrium. We shall be satisfied in this chapter with formal derivations 
(for the rigorous treatment of uniformly hyperbolic dynamical systems, see Chapter 5). The 
rest of this chapter is based in good part on Gallavotti and Ruelle [46] . 

2. Linear response close to equilibrium. 

Given an SRB state p for (/*), we have introduced in Chapter 4 the response function 
a I— > K a such that K a = for a < and 

(a)$ = J p{dx){{T x{ _ a) r)x{r a x))-v x § 

= J p(dy)X(y) ■ V y (d> o n for a > (3) 

With this definition we have shown that the linear change 5 t p of the SRB state p in response 
to a change 5 t F in the force F acting on the system is given by 

&tp= J do-K a 5t- a F (4) 

in the continuous time case (t e R). In the discrete time case (t G Z) we have with the 
same definition of K a 

5 tP =J2^t-afof- 1 ) (5) 
a 

where S t f is the infinitesimal change of / at time t (and S t f o f~ x is thus an infinitesimal 
vector field on M). 

We shall denote by AivV the divergence of a vector field V on M with respect to the 
volume element p(dx). By definition, the Lie derivative of this volume element satisfies 

L v p = (divV)p 

We may thus rewrite (3) as 

faX)* = J p(dy)X(y) ■ V y (d> o r ) = - J p(dy)drvX(y)Mry) (6) 

for a > 0. 
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For the susceptibility we obtain in the continuous time case (w e R): 



(k w X)$ = e^da J p{dy)X{y) ■ V„($ o f a ) 

= - e^da J p(dy)divX(y)Mry) 
and in the discrete time case (oo G R(mod 2n)) 

oo „ 

fox)* = ^ / p(dy)X(y) • V y (d> o r) 



a=0 



oo „ 

/ p(dy)divX(y)MFy) 



(7 = 

In particular for time independent 5F or 5f we have respectively 

coo 



(6p)$ = (k 5F)® = - da p(dy)dw5F(y)Mf a y) 



o 



or 



(6p)* = (k (Sf o = - E / p(^)( div W ° /- 1 ))(2/)-$(/ fJ 2/) 

cr=0 

3. Entropy production. 

Suppose that p is an equilibrium state corresponding to the vector field F (or the 
diffeomorphism /). For a small time dependent change 5 t F of F (or S t f of /) we shall be 
able to determine the entropy production to second order. 

Remember that in a time independent situation we can (as noted in Section 1.9) use 
any smooth volume element to compute the divergence (or the Jacobian) occuring in the 
definition of the entropy production. To estimate a time dependent entropy production 
we have however to make a definite choice of volume element. We take here p(dx), which 
is the natural choice. In the continuous time case we have then 

divF = 

[Indeed, by definition, p(dx)divF(x) = Lpp(dx) = because p is invariant under the time 
evolution defined by F] . In the discrete time case we have 

J = 1 

[because / preserves p(dx)]. We expand the SRB state corresponding to F+S t F (or f+8 t f) 
to second order as p + S t p + 5^ p, and compute the corresponding entropy production. 
The term corresponding to 8^ p does not contribute to the result, which is as follows. 
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Proposition (time dependent entropy production). 

In the continuous time case, writing X t = S t F, the entropy production is to second 
order 

et = L da I p( d y"> dwXt -M- dwXt (f a y} ( ? ) 

In the discrete time case, writing X t = 5 t f o f~ x , the entropy production is to second order 

e t = \j p(dy)(divX t (y)) 2 + J2J p(dy)&vX t - a (y).divX t (ry) (8) 

In the continuous time case we have to second order 

e t = - J (p(dx) + S t p(dx) + 4 2) /9(^))(divF(x) + div5 t F(x)) 

Since divF = 0, the 5^ p-term does not contribute, and since the p-integral of a divergence 
vanishes we find 

e t = - J S t p{dx)div8 t F(x) = -(8 t p)(divX t ) 

Using the linear response (=first order) formulae (4), (6) gives then (7). 

In the discrete time case, for the entropy production between time t — 1 and time t, 
we have 

e* = - J (p(dy) + 5 t -ip{dy) + &?\p{dy))\ogJ f + 6tf (y) 

where 

J f+5 t f(y) = Jidentit y +5 t fof-i(fy) = 1 + higher order 

so that the 5( 2 )p-term vanishes. Let us now suppose that the vector field X t = 5 t f o / _1 
has support covered by a local chart on which p is proportional to Lebesgue measure: 
p(dx) = cdx. Using the coordinates of this chart, and denoting by A the matrix (diXj), 
with X = X t , we have to second order 



p(dy) log Jf+S t f(y) = J P{dx) log J id entity+X t (a^) 

= cj e£rlogdet(l + A) = c J dxtrlog(l + A) 
= c J dxtx{A- l -A 2 )=c J dxiY^diX^-^idiXjix^djX^x))) 

i ij 

Using integration by part, this is 

= ~\f dxiJ^Wx)) 2 = - 1 - I p(dy)(^X(y)f 
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For general X we may write X = J2X a and assume that every |X Q | + \Xp\ has support 
as above; from this one obtains that the formula 



p(dy)\ogJf+ St f(y) = ~ J p{dy)(divX t (y)f 



holds generally (to second order). 
Also to second order we have 



J dt-ip(dy) log J f+St f(y) = J 5 t -ip(dy) log Ji d entit y +x t (fy) 

and therefore 

e t = \\ p(dy)(divX t (y)) 2 - [ 5 t _ lP (dy)divX t (fy) 



2 _ 

Using the linear response (=first order) formulae (5), (6) gives finally 

e t = \j p(dy)(divX t (y)) 2 + J2j p(dy)divX t _ a (y).divX t (r y ) 

which proves (8). [j 

Corollary (time independent entropy production). 
If X does not depend on t, 

1 f°° f 1 f°° 

€= 2J da J P( dx ) diyX (y)- diyX (f a y) = 2 J da Pdwx,di v x((T) 

in the continuous time case, and 

Yl / P( dx ) divX (y)- divX (f a y) = g Yl Pdivx,divx((7) 

cr= — oo cr= — oo 

in the discrete time case. 

[Remember that p AA (a) = f p(dy)A(y)A(f a y) - (f p(dy)A(y)) 2 J p{dy)divX{y) = 0]. 
Remark (frequency dependent entropy production). 

In the continuous time case, taking X t (y) = \pl sin {ut — 6)X(y) and averaging over 
t we obtain ^ 

(e)(u) = J^ cosuada J p(dy)divX(y).divX(f CT y) 
1 f 00 

= e tuja dapdi v x,divx(v) 

Z J -oo 
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= ^PdivX,divx(w) = 7T/5divX,divX (^) > 



[The right-hand side is proportional to the spectral density associated with the signal 
divX o / and is thus > 0, see Section 4.9]. 

4. Forces and observables close to equilibrium. 

In the last two sections we have studied the linear response and entropy production 
corresponding to a force X (= 5F or 5 f o / pushing a system away from equilibrium. A 
striking feature of the formulae that we have obtained is that X occurs only in the form* 
divX. Nonequilibrium statistical mechanics close to equilibrium is very special because of 
this feature, which does not persist far from equilibrium. 

Let us call observable a function $ : M — > R. We further assume that $ is smooth 
and that p($) = 0. Thus, if X is any smooth vector field, divX is an observable, and 
every observable is of the form divX. [More precisely, if M is connected, every function 
with a vanishing integral can be written as a divergence]. From a physical viewpoint, it is 
essential to be able to identify the observable divX correponding to a force X . Remember 
that we typically deal with macroscopic systems, for which the mathematical definition of 
the divergence has no direct operational meaning. 

To identify divX as an observable we shall use the Hilbert space 



The observables $ belong to Ti and the scalar product (^J, $) is = p*$(0) (value 

at of the correlation function p^$(t). We have in principle an operational definition of 
the scalar product 



Indeed, by (4), (6) this is the response at time 0+ of a time dependent force Xt = 5(t).X 
where d(t) is the Dirac delta in the continuous time case (the discrete time case is similar). 
The force X t is a kick at time 0, and (— divX, $) is the instantaneous response of the 
system**. One can also say that the kick X t replaces p(dx) by (1 — divX) p(dx) at time 
0+, where the new term — div X p(dx) may be interpreted as a fluctuation. [We do not 
wish to pursue the discussion of fluctuations on the basis of Markovian and Gaussian 
assumptions as done for instance in [20]]. 

* Remember that divX is the divergence of X with respect to the volume element p(dx) 
corresponding to the equilibrium state p. 

** In the continuous time case, one could also switch on a constant force X at time 
and obtain (— divX, $) as derivative of the expectation of $ at t = 0+. In the discrete 
time case, (— divX, $) is the jump in the expectation value of $ when / is replaced by 
(identity + X)of. 



W={$GL 2 (M,p):/)($) = 0} 
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Notice that we can also obtain the scalar product (divX, divF) by polarization of 



(divX,divX) = p((divX) 2 ) = i J (e)(uj)doj 



The formulae written above permit in principle to identify divX as an observable. 
From now on we accept without further discussion this identification, which is important 
in particular for the Onsager reciprocity relations and the fluctuation-dissipation theorem. 

5. Thermodynamic forces and conjugate fluxes. 

Let us assume that the dynamics (/*) of our system depends (smoothly) on parameters 
E a , as is natural in physical applications. We also assume that when the parameters E a 
take the value equilibrium is achieved, corresponding to a vector field F (continuous time 
case) or a diffeomorphism / (discrete time case). Close to equilibrium F is replaced by 
F + 5F, or / by / + 5f, and we take as earlier X = 5F or X = 8f o We may thus 
write, to first order 

X = Y,V a E a 

a 

{i.e., SF = ^2 a V a 5E a or Sf o / _1 = ^2 a V a SE a , since E a = SE a ). We assume here 
that the E a are time independent; the time dependent case is discussed in Section 6. The 
thermodynamic forces which occur in the formalism of nonequilibrium thermodynamics 
(see [74], [20]) will now be identified with the E a . 

In the continuous time case we define, following Gallavotti [37], [38], the thermody- 
namic flux conjugate to E a as 

T ( 9 \ 

where 

a F+SF = -div(F + SF) = -divX 
is the local entropy production. To first order in X the flux is (using (4), (6)) 

= 5p(-divV a ) = J da J p(dy) divX(y).divV a (ry) = ^ LapEp (9) 

where 



Lap = J q da J p(dy) d\Wp(y).d\WaU a y) 



In the discrete time case we define the thermodynamic flux conjugate to E a as the 
following first order expression in X, analogous to (9), 

x a = l p (divxdivvg - d P ((divv a ) o /) 
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2 

where 



OO 

\ [ p(dy) d\vX{y).d\W a {y) + ]T / p(dy) divX (y) .divV a (f a v) = E L «^ ( 10 ) 

1 J a=l J 



L ^ = \j p(dy)diyV (3 (y).diyV a (y) + f^J p{dy)diWp{y) .dxvV a {r y) 

The definitions (9) and (10) are such that the entropy production ( given by (7) or 
(8)) is 

a a/3 

6. An operator formulation. 

Let us assume that there is a (sufficiently large) Banach space B of continuous func- 
tions $ : M — > R such that p(G>) = and 



/: 



da\ P (*.(<s>or))\<c\\<f>M\nB (ii) 



or 



da\p(*.(*or))\ <C||$|| b ||*||b (12) 

cr = — oo 

if $, \1> G B. In particular, (12) holds if / is an Anosov diffeomorphism and B is the space 
of a-Holder continuous functions for some a > 0, with vanishing p-integral. 

We want to take advantage of the fact the formulae for J a only involve the divergence 
of X and V a to give an operator formulation of the relation between thermodynamic forces, 
fluxes, and entropy production, with time dependent Xt = 5tF or 5tf ° / • Assuming 
that the divergence of X t is in B, we define X t G B and X t G £>* (the dual of B) as follows* 

X{t) = -d\vX t 

and 

X(t) — Stp (continuous time) 

or 

1 

X(t) = — — {divXt).p+ f*5t-ip (discrete time) 

[In the time independent situation of Section 5 we have thus X a = (X(0), — divV a )]. The 
entropy production (7) or (8) may now be written 

e t = (X(t),X(t)) 



* We give here a corrected version of (8), (8') of [46]. 
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Let us introduce a linear operator K(a) : B — > £>* such that 

{0 for a < 
§p(tf$) for a = 
^.(Sof )) /or a>0 

[This is a slight modification of the definition of K a ; K(a) is a response function adapted 
to the situation close to equilibrium. Note that the value at a = is relevant only in the 
discrete time case]. We have now 

l(t) = J daK(a)X(t-a) 

or 



X{t) = Y J K{a)X{t-a) 



in the discrete and continuous time cases respectively. Note that causality is respected in 
the sense that X(t) depends only on X(t — a) for a > 0. Let L{u>) : B — * B* be the real 
part of the Fourier transform of K(a): 



or 



poo 

L(u) = da cos ua K(a) 
Jo 



= ^ coscua K(a) 



cr = 



[The convergence of the right-hand sides is ensured by (11) or (12)]. In the time indepen- 
dent case (writing X(t) = X, X(t) = X) we have thus 

1 = L(0)X 

More generally, let T{u>) be defined (in the continuous time case) by 

poo 

(IM,$) = ( / dap((-divX t _ a ).<f> t or)) t 
Jo 

where X t = v / 2sin(o;t - 9)X, $ t = V2sin(ut - #)$, and (. . .) t is average over t; taking 
X = — divX we have 

X{u) = L(uj)X 

The discrete time case is similar. 

7. Onsager reciprocity and the fluctuation-dissipation theorem. 

We may interpret the Fourier transform K(uj) of the response function K(a) as the 
susceptibility. We have in the continuous time case 



poo 

Jo 
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and in the discrete time case 

1 oo 

CT = 1 

In this section we shall make use of reversibility. We may thus define an operator 
e : B — > i3 by e$ = $ o z, and we have e 2 = 1. Also B = B + + B- where £> + = {$ + e$ : 
$ G i3_ = {$ - e$ : $ G B}. Note that if we write e($) = ±1 when $ G #± we have 
e $ = e($)$. 

Proposition. 

(a) (K(u)^ : $) extends analytically to Imw > (by causality). 

(b) (K(uj)-V, $) = e($)e(#)(K(u;)$, *) if $ G B±, * G £± (by reversibility). 

(c) The real part L(u) of K{uj) satisfies 

(L(u)V,<$>) = p^{u) (13) 
e* L(u)e = L(u) (14) 
Part (b) and equation (14) result from 

p(*.$ o D = (i*p)(M o D = p((vp o z).($ o i o /-)) = ^.(e*) o n 

The other assertions are immediate. [] 
Remark. 

By (13), the quadratic form on B defined by $ i— > (L(cu)$, $) is > because /5$$ > 
as discussed in Section 4.9. 

Reciprocity relations (Onsager [68]). 

Taking uj = and writing L(0) = L, so that X — LX in the time independent case, 
we obtain from (14) the usual reciprocity relations 

(L^,$) = e($)e(#)(L$,#) 

or, with e a = e($ Q ), ep = ei&p), and the notation of Section 5 

La/3 = tottpLpa 



Fluctuation-dissipation theorem (Green [50], Kubo [57]). 
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The formula (13) gives 

±{(K(u)(-divX), $) + (K(-u)(-divX), d>)] = p { - divX) ,^) 

i.e., the real part of the susceptibility corresponding to the force X is given by the Fourier 
transform of a time correlation function for the equilibrium state p. This is a form of 
the fluctuation-dissipation theorem. The left-hand side corresponds to dissipation (the 
nonequilibrium process associated with X) and the right-hand side to fluctuation (in the 
equilibrium state p). The physical use of the above formula requires that we can identify 
— divX as an observable (see Section 4). We have discussed in Section 4.10 what remains 
of the fluctuation-dissipation theorem far from equilibrium. 

8. Higher order corrections, and noisy systems. 

Consider the continuous time dynamical system (/|) associated with the vector field 
F + XX on M. [It is convenient here to make the powers of A explicitly visible, while in 
earlier sections we took A = 1. We take F + XX to be of degree 1 in A for simplicity]. 
Given an SRB state p for (/q) and assuming that one can perturb p to an SRB state p\ 
for small |A|, we obtain (see [87]) a formal expansion of p\ in powers of A, given by 

°° poo poo 

p A ($) = VA r / dr x ... dT r p (PQ Tr P...PQ Tl $) (15) 

where $ : M — > R is smooth and independent of A, and the operators P, Q l are defined 
on smooth functions by 

P^> = X ■ V* 

Qty = ^of t 

Formally, we can also write 

poo 

p x ($) = p ((l-X / drPQ^y 1 ®) 
Jo 

The simplicity of these formulae is deceptive. In particular, if one decomposes X into 
stable, unstable, and neutral components to study the convergence of the multiple integrals 
in (15), one gets complicated expression which appear unusable in practice. If one assumes 
that po is an equilibrium state one recovers to leading order the results of Sections 2 and 3 
for Sp and e, but the higher order contibutions can no longer be expressed solely in terms 
of divX. 

It is interesting at this point to discuss systems with a little bit of noise added. We start 
with a probability space (O, P) and a continuous one-parameter group (#*) of transforma- 
tions of O with respect to which P is ergodic. A one-parameter group of transformations 
on Q x M is given such that the time t map is (u>, x) i— > (6**a;, f^x). One could here again 
define SRB states (see the discussion of the discrete time case in Section 4.2). We shall 
however take another approach, forget about hyperbolicity, and discuss the Markov case 
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which is realized by stochastic differential equations (see for instance Arnold [2]). Let 
(ft\) ^ e the family of diffeomorphisms obtained by integrating the stochastic differential 
equation (Langevin equation) 



dx 
~dt 



F(6 t u,x) + \X(x) 



and define 



= X ■ Vtf 



If p is a probability measure with smooth density on M we define a stationary measure 
px as a weak limit of J P(du)(f^ x )* p when T — > oo. Using the Markov property we find 
then formally 



and p T is the kernel of the diffusion process associated with the Langevin equation. 

We must refer to [87] for a more detailed study of this problem, but we note the sim- 
ilarity of the expression obtained for the stationary measure p\ and the earlier expression 
obtained for the SRB measure p\ in the deterministic case. In fact, in the zero noise limit, 
the coefficients of the A expansion for the stationary measure formally reproduce the coef- 
ficients for the SRB measure. This is in accordance with various results on the stability of 
SRB states under small random perturbations for uniformly hyperbolic systems (see Kifer 
[55], [56]). Since the roundoff errors in numerical simulations may be viewed as a small 
random noise, it is satisfactory that they respect SRB states. 

We again refer to [87] for the study of higher order corrections for discrete time 
dynamical systems. 




where 
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Chapter 6. ISOKINETIC TIME EVOLUTIONS 
AND THE DETTMANN-MORRISS PAIRING THEOREM. 



1. Isokinetic models. 

Consider the evolution equation 



— ( p \ = ( ^(q) — «p 
dt 1 q J V p/ m 



(i) 



where p, q G R^ and a = p • £/p • p. This is called an isokinetic time evolution because 
it conserves the kinetic energy p • p/2m, and we may thus restrict (1) to the manifold 
S x R^ where 

S = {p g K N : p • p/2m = K} 

We may replace the configuration space R^ by an open subset M of R^ or of the 
iV-dimensional torus T^, with a piecewise smooth boundary dM. The time evolution (1) 
in S x M is then complemented by elastic reflection at the boundary dM (see Section 5). 
Such a model is discussed in [16], [17] (where £ is taken to be the sum of an electric force 
and a magnetic force) . A system of n hard balls in R d can also be treated in this manner. 
Here we take N = nd and, if the balls have different masses mjt, positions x k and velocities 
v k , let 

( m k \l/2 / m \l/2 

m m k 

(with these variables p, q, the shocks between the balls correspond to elastic reflections at 
the boundary dM). 

A generalization of the isokinetic evolution (1) consists in replacing the configuration 
space R^ by a manifold M with a Riemann metric (see Section 10). In the absence of 
force we obtain thus a geodesic flow on M. 

Suppose that the force £ in the isokinetic evolution equation (1) is locally a gradient 
(this is the case for an electric, but not for a magnetic force). Then Dettmann and Morriss 
[22] have proved a pairing property of Lyapunov exponents (see Theorem 7) . Such a pairing 
has been observed in the numerical study of various models (see [21], [61]), and its rigorous 
justification (in [22], [23], [24], [18], [96]) is the main reason for the interest enjoyed by 
isokinetic models. 

2. Entropy production. 

The divergence divF of the right-hand side of (1) is readily computed to be —(N — l)a 
(this does not require that £ be locally a gradient). The entropy production e = e(p), which 
is equal to minus the expectation value of divF, takes thus the simple form 

N - 1 

e{p) = {N-l)p{a) = -^w 
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where w = p(p-£/m) is the work performed on the system (or dissipated) by unit time, and 
N — 1 may be interpreted as the number of degrees of freedom. In equilibrium statistical 
mechanics the "equipartition of energy" would give K/(N — 1) = T/2. (We take here 
Boltzmann's constant k = 1, as we have done for the definition of entropy production). 
Therefore 

w 

6= T 

which fits the formula 

T 

where S is the entropy, but notice that here we may be far from equilibrium. 
3. The variations 5p, 5q. 

We suppose p ^ and consider a solution X : t t— > (p, q) e S x R^ of (1), and an 
infinitesimally close solution Ja : t i— > (p + Ap, q + Aq) also in S x R^. Define Sp = Ap 
and Sq = Aq — (p • Ap/p • p)p, then p ■ dp = p • Sq = 0, and X$ : t i— > (p + 5p, q + Sq) 
is (neglecting higher order terms) a reparametrization of It will be useful to consider 
Xg when we discuss Lyapunov exponents for the flow defined by (1). 

Lemma. 

Let (e^)^ 1 be an orthonormal basis of the subspace of H N orthogonal to p, so that 
dei/dt = -(£ • e;/p • p)p. Writing Sp = Y^iLi 1 $Pi e i, = Y.iLi 1 % e i, w ? ^ave 



—i , , , (2) 

dt \5qi J \(l/m)8ij J \5qj J 

Note that the conditions on the dei/dt ensure that (ej) remains an orthonormal basis 
of the subspace of Tt N orthogonal to p (they define parallel transport of (e$) along X). 
We have 

- Ap = (Aq • d)£ - aAp - (Aa)p 



dt 



d . 1 . 
-Aq= -Ap 

dt m 



We may replace Ap, Aq, by Sp, 5q, in the right-hand sides (Aq is initially chosen orthog- 
onal to p), then take a projection parallel to p to obtain (2). [] 

4. Proposition. (Dettmann-Morriss [22]). 

If Sx(t) is the vector with components Spi, 5qi, we may write 

8x(t) = L(t)8x(0) 

where L(t) is a (2N - 2) x (2iV - 2) matrix and L(0) = 1. Also let 

Si 



J 



-Stj o 
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Then if £ is locally a gradient, we have 

L(t)*JL(t) = exp(- I a{X{s))ds)J 
Jo 

where L* is the transpose of L, and a(x) = p ■ £/p • p = p • £/2K when x = (p, q). 
In view of (2) we have 

±L { t) = TL(t) 

where 

( ad,, diij \ 

so that 

T*J + JT= ( ° _ " a<y « ^ 
If £ is locally a gradient this gives 

T*J+JT = -aJ 

hence 

^[I(t)Vl(t)] = -a[L(t)VL(t)] 
and since L(0)*JL(0) = J, we obtain 

L(t)*JL(t) = /iJ (3) 

with n = exp(— J Q * a(X(s))ds). [] 

The formula (3) is expressed by saying that L(t) is conformally symplectic. 
5. Proposition. (Elastic reflection). 

Elastic reflection at an obstacle corresponds to a linear transformation 5X — > 5 R X = 
R5X where R is symplectic (R*JR = J). 

This result, in one form or another, is very well known. For completeness we give a 
proof. As in the previous proposition, 5X e H 2N ~ 2 is a vector with components Spi, dqi 
such that dp = Yl,i&Pi e i-> = Y^i&Q^ii an ^ p • 5p = p ■ 5q = 0. Similarly for 5 R X: 
$ R P = S R Pi e R i, 5 R ci = £\ S R qi e R i. 

Notice that 5 R q can be computed (to first order) as if the obstacle were planar, with 
normal vector n, so that 

S R q = (5q-2(n- 5q)n 

Taking 

e R i = ei - 2(n • e;)n 
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we have thus 5 R qi = 5qi for i = 1, . . . , N — 1. 

To obtain 5 R p, we may assume that the surface S of the obstacle is given in by 
xn = <&(xi, . . . , xjv-i) where $ and its first order derivatives vanish at (0, . . . , 0). Neglect- 
ing higher order, we see that the normal n + 5n to E at (xi, . . . , xn-i, $(^i, • • • ,xn-i)) 
is (-<9i<E>, . . . , -9at_i$, 1). We have 

p R + S R p = p + dp - 2[(n + <5n) • (p + <5p)](n + <m) 

so that 

5 R p = 5p- 2(n • p)<5n - 2[(n • <5p) + (p • <5n)]n 

= 5p- 2(n • p)5n + 2(p • <5n)n - 2[n • (5p - 2(n • p)<m + 2(p • <5n)n)]n 

With our choice of (e R i), this means that the components of 5 R p with respect to (e R i) 
are 

5 R pi = dpi - (2(n • p)5n - 2(p • <m)n) • e< 
or, writing 5p* = £V & R Pi e i, we have 

5p* = <5p - 2((n • p)<m - 2(p • <5n)n) 
= dp — 2(n ■ p)Q5n 

where Q = -P -1 is the inverse of the orthogonal projection P of the hyperplane {u : n - u = 
0}cR w to {u:p-u = 0}cR w . It is readily seen that 

8n = -AQdci 

where A is the matrix if second order derivatives of $ at 0. Therefore 

Sp* = 5p + 2(n • p)QAQSq 

The matrix R such that 5 R X = R5X is thus 

R= ^l 2(p-n)QAQ^j 

Since the matrix A is symmetric, QAQ is symmetric, and R*JR = J, proving that R is 
symplectic. [] 

6. Lemma (Spectrum of a conformally symplectic matrix). 

If L*JL = fiJ, then the spectrum of L*L (including multiplicities) is invariant under 
A I— > /U 2 A _1 . Tie eigenspaces for the eigenvalues A and // 2 A _1 are interchanged by J. 

By assumption 

L* = /j,JL~ 1 J~ 1 
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Using also the fact that J is unitary we obtain 

L*L = p 2 JL~ 1 L*~ 1 J~ 1 = p 2 J(L*L)~ 1 J~ 1 

This unitary equivalence proves the lemma since J -1 = —J maps the eigenspace for the 
eigenvalue A to the eigenspace for the eigenvalue p 2 X~ 1 . [] 

7. Theorem. (Pairing of Lyapunov exponents). 

Let M be an open subset of~R N or of the torus T N , with a piecewise smooth boundary 
dM . Consider the time evolution onSxM defined on SxM by (1) with £ locally a gradient, 
and at dM by elastic reflection. Assume that p is a probability measure on S x M such 
that f f x is well defined p(dx)-a.e., and that p is ergodic under the time evolution (/*). 
Then, the 2N — 1 Lyapunov exponents associated with (/*) and p may be labelled Aj for 
i = —N + 1, . . . , N — 1 in such a way that 

A = , A; + A_i = -p(ct) for (4) 

(with p(a) = f p(dx) a(x)). 

There is a zero Lyapunov exponent associated with the direction of flow, the other 
2N — 2 are obtained as logarithms of the eigenvalues of 

A= lim(L(t)*L(t)) 1 / 2 * 

for p-a.e. x, and L(t) defined by 8x(t) = L(t)8x(0). Writing L(t) = L, we have by 
Propositions 4 and 5 

L*JL = pj 
with p = exp(— J * a(x(s))ds). Therefore 

L = -pJL*~ x J 

L*L = p?J- 1 (L*L)- 1 J 

so that 

A = (exp-p(a))J- 1 A" 1 J 

Using Lemma 6 we see that the eigenvalues of A come by pairs (u,v) with uv = e~ p ^; 
the theorem then follows readily. |j 

Remark. 

With the above notation there are subspaces V-n+i C ... C V_i C V\ C ... C 
Vat-i = R 2N ~ 2 such that 

lim - log ||L(t)w|| = Xi for u G VAV- 

t^oo t 
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where we have written V- = Vj-± if j 7^ 1 and V[ = V-\. The spaces Vj are such that for 
j > 1, V-j is the skew-orthogonal J^V-) 1 - ofVj. 

If Ui is the eigenspace of A corresponding to the eigenvalue exp(A^) we have Vj = 
®i<jUi and, since U-j = JUj, 



r\±- 



V-j = ® % >jU-i = J ®i>j Ui = J(Vj) 

which establishes the asserted skew-orthogonality. [] 

[Basically, Theorem 7 is due to Dettmann and Morriss [22], but the inclusion of shocks 
and the above remark on skew-orthogonality are due to Wojtkowski and Liverani [96]]. 

8. Time dependent forces. 

Let P be an ergodic measure for the flow (6>*) on O, and let £ w be a vector field on 
M depending on w 6 11. Writing also ct w = p • £ u /p ■ p, we obtain a time evolution 
X ( Q ) fi x (°) = x (t) = (P> q) b y solving 

^_ fp\ = ( £o*u (q) -®etu;P 
dt\(l) \ p/m 

on S x M, with elastic collisions at dM. Clearly f^ +s = /|t w /w, so that (a;, t) h-> f*(a;, t) = 
{O 1 ^, fljX) defines a flow (f*) on O x S x M. The formula for the entropy production is 
here 

e(p) = (N — 1) J p(dudx)a u (x) 

If p is an (f *)-ergodic measure (with projection P on O) we can define Lyapunov exponents 
Xi such that the exp are the eigenvalues of 

A., x = hm((T x ^)*^)) 1/2t 

for p-almost all (u>,x). If ^ is locally a gradient (for P-almost all u), then the Aj again 
satisfy the pairing property (use the same arguments as in the time independent case). 
Here 

Xi + X-i = - J p{dudx)a UJ {x) 

9. Generalizations: conformally symplectic structures. 

Let M. be a 2iV-dimensional manifold. Wojtkowski and Liverani [96] define a confor- 
mally symplectic structure on M to be a differentiable 2-form G which is non-degenerate 
and such that 

dQ = 7 A e 

for some closed 1-form 7. Locally, we have thus 7 = dU and d(e~ u <3>) = 0. 

Given a function H : M. — > R (Hamiltonian), we define a vector field F = VqH by 

Q(-,F) = dH(.) (5) 
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Since e u F is locally a Hamiltonian vector field with respect to the symplectic form e ^O, 
the conformally Hamiltonian flow (/*) defined by F conserves the energy H. 

Let M K = {x G M : = and 6 K be the restriction of 6 to T m M K - If 

it, u G T x A^k , we write 

(^e*)(«,«) = |e(T x /VT x A)| t=0 

Theorem. (Wojtkowski-Liverani [96]) 

For a conformally Hamiltonian How restricted to the level set Mr, we have 

jQ K =l{F)Q K 

where F is defined by (5). 
On Mk we have indeed 

d{e~ u {H - K)){-) = e~ u dH{-) = e" c/ e(-, F) 

i.e., F coincides with the Hamiltonian vector field given by the Hamiltonian e~ u (H — K) 
with respect to the symplectic form e~ u Q. Since the symplectic form is preserved by the 
Hamiltonian flow we have 

o = i(e- u e K ) = -e-^e, + e-^e, 

dt dt dt 

where 

^ = dU(F) = 7 (F) 

and the theorem follows. [] 

Since Q(T X M k, F(x)) = 0, 0^ defines a bilinear form on T x Mk /R~F(x), and we 
may choose a basis of this (2N — 2) -dimensional space so that 6^ is defined by the matrix 
J. Therefore the theorem implies that T x /*(mod R.F) is given by a matrix L(t) such that 

L(t)*JL(t) = J.exp f 1 {F{f t x))dt (6) 
Jo 

We shall now show how to recover Theorem 7 from (6) in a somewhat more general 
situation*. 

10. Isokinetic dynamics on a Riemann manifold. 



* Another application treated by Wojtkowski and Liverani [96] is to Nose-Hoover dy- 
namics (first discussed by Dettmann and Morriss [24]). 
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Let M be an iV-dimensional Riemann manifold, and £ = ^2iLi£i(<l)dq l a closed 1- 
form (i.e., is locally a gradient). We let M. = T*M be the cotangent bundle of M; its 
points are of the form (p, q) with q G M, p G T*M. We take as Hamiltonian the kinetic 
energy 

where {g 1 ^) is the inverse of the matrix (gij) defining the metric. Write also 



K 



= ^^Pidq 1 , u = dn = dpi A dq l 



e = u + ^ AK 

Then 

dO = A dn = A e 

2K^ 2K^ 

i.e., © is conformally symplectic. 

With respect to the basis (dp i7 dq l ), and dH are given by 

Sij \ / g ij pj/m 

Sij [ZiPj - ZjPi]/2K J ' \W k /dq i )p jPk /2m 

so that the conformally Hamiltonian evolution equation is 

±(Pi\ =F= ( [(9 ik PjPk)^i - (g j %Pk)Pi]/2mK - (dgi k /dq*)p jPk /2m 
dt\q*J V 9 ij Pj/m 

with the usual summation convention. On the level set M.k this reduces to the expected 
analogue of (1): 

^_(Vi\ = (ii- {g ik ijPk/^rnK)pi - (dg jk /dq l )p jPk /2m\ 
dt\q l J V g tj Pj/m J 

The pairing of Lyapunov exponents follows thus from (6) and Lemma 6: 

A = , A; + X-i = -p(a) 

where a = —7(F) = g l ^iPj/2mK. If M has a piecewise smooth boundary and the dy- 
namics is complemented by elastic shocks at dM, the above considerations remain correct 
because elastic reflection is symplectic. 

[Pairing for isokinetic dynamics on a Riemann manifold has been established by 
Choquard [18]; the above treatment follows Wojtkowski and Liverani [96]]. 

11. Homological formula for p(a). 
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If £ is locally a gradient, i.e., if the 1-form £ = £ ■ <iq is closed, it defines a cohomology 
class [£] G if 1 (M). Let C(x,t) denote the curve {q(s) : < s < t} associated with the 
initial condition x = (p(0),q(0)), then the limit 




C(x,t) 



exists for all continuous 1-forms uj, and is p(dx)-a,.e. constant (by ergodicity of p). This 
is a linear functional of uj, which vanishes on exact forms (i.e., when uj = d<&) and defines 
thus a homology class [p] e H\(M). We may thus write 



12. Thermodynamic limit. 

In equilibrium statistical mechanics "large systems" exhibit thermodynamic behavior. 
This means in particular that when the total energy is fixed, the kinetic energy has only 
small fluctuations around some average value, and vice versa. It has thus been suggested by 
G. Gallavotti (see [30], [40] ,[41], [42]) that for a large system the isokinetic time evolution 
defined by (1) would be equivalent to a (more physical) time evolution with fixed total 
energy. 






([£]>]) 
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Chapter 7. LOOSE ENDS AND OUTLOOK. 



1. Frameworks for the study of equilibrium. 

The view of nonequilibrium statistical mechanics presented in previous chapters was 
based on chaotic microscopic dynamics, and centered on nonequilibrium steady states. We 
have given priority to technical results (rigorous and less rigorous) over general ideology. 
In particular we have restricted ourselves to classical systems, having little to say about 
the quantum case. At this point however it is desirable to open up the discussion to more 
general possibilities. 

We have chosen to analyze stationary nonequilibrium situations, rather than study 
the approach to equilibrium a la Boltzmann. While maintaining this general philosophy, 
we shall see that the restriction to steady states is artificial and must be somewhat relaxed 
(Section 2). 

To achieve both nonequilibrium and some kind of stationarity, we are led to discussing 
either infinite systems or nonhamiltonian forces. We shall now briefly review these two 
cases, disregarding other approaches, in particular the important escape rate method of 
Gaspard, Nicolis and coworkers [48], [47]. 

Infinite systems. 

It is customary to think of a finite system E interacting with infinite heat reservoirs. 
The reservoirs are often described as free fields at different temperatures. Since the in- 
teraction of the free fields with E is localized, those degrees of freedom which have been 
influenced by E conveniently move away to infinity and can be forgotten. The case of a 
single reservoir is used to study approach to equilibrium (see Jaksic and Pillet [53]). To 
achieve a nonequilibrium steady state, one needs at least two reservoirs: 



The study of this situation is still in an early stage (see Rieder et al. [77], Lebowitz 
and Spohn [94], [63], Eckmann et al. [29], [30]). Note that in a macroscopic description of 
a realistic reservoir in a steady state, the temperature T would satisfy the heat equation, 
hence (with constant heat conductivity) AT = 0. A limiting value T\ or T-i at infinity is 
thus possible only if the space dimension is > 3. This shows the importance of diffusive 
effects in nonequilibrium, and how an oversimplified model can easily give incorrect results. 
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In the infinite system situation now discussed (with Hamiltonian forces) there is no 
compelling reason to consider classical rather than quantum systems. Note also that dissi- 
pation here occurs because information moves away at infinity in the reservoirs and is lost. 
The study is much simplified when free fields are used as reservoirs because the degrees 
of freedom in the reservoirs are independent (no rescattering). Incidentally, the absence of 
rescattering is also the key reason for the success of Lanford's analysis of the Boltzmann 
equation in the Grad limit [60]. Studying the loss of information at infinity in the pres- 
ence of rescattering seems to be a formidable problem. Since a realistic understanding of 
nonequilibrium in infinite systems depends on the solution of this problem, it is natural to 
"cheat" and try to understand nonequilibrium in finite systems. 

Finite systems with nonhamiltonian forces. 

These are the systems studied in previous chapters. A typical example is the following 
system E: 




E consists of classical particles confined to a solid torus, and having ordinary (Hamil- 
tonian) interactions with each other and with scatterers. Furthermore the particles are 
charged and pushed around the torus by an electric field S associated locally with with 
a potential U. In fact the electric potential is multivalued, and changed by an additive 
constant after one turn around the solid torus. Finally, the system is thermostatted by 
fixing the total kinetic energy of the particles. The multivalued electric potential is only 
mildly objectionable from the point of view of physics, and so is the isokinetic thermostat 
(which could be replaced by another similar thermostat). Altogether one has a reasonably 
realistic physical system with the great virtue that, as we have seen, it allows a nontrivial 
theoretical treatment. 

One may well argue that in statistical mechanics (nonequilibrium as well as equilib- 
rium) the infinite system limit (thermodynamic limit) must be be taken. We shall look 
into this problem in Section 3. But note here that in equilibrium statistical mechanics it 
is advantageous to put the study of ensembles describing finite systems at the beginning, 
and leave for later the discussion of the thermodynamic limit, or the analysis of actually 
infinite systems. A similar strategy seems reasonable also for nonequilibrium. 

Unfortunately, the nice results obtained for finite classical systems with nonhamilto- 
nian forces seem difficult to extend to the quantum situation. Indeed, the Schrodinger 
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equation requires a true potential function, not a multivalued one*. One may also remark 
that if SRB states can forever produce entropy, it is because their own entropy is — oo. A 
quantum analog of this situation seems hard to achieve, because quantum entropies are 
positive. 

2. Nonsteady nonequilibrium states. 

We have taken as steady nonequilibrium states the limits 



where m is a probability measure absolutely continuous with respect to the Riemann 
volume element dx on the phase space M. But we have also asked that p be ergodic or, 
better, sufficiently rapidly mixing (this is needed in applications like the linear response 
formula). Of course one can always decompose an invariant state p into ergodic components 
but these need not be of the form (1), or they may be of the form (1) but not mixing**. 

It is not clear how to address this problem, which appears related to the question 
of macroscopic oscillations of nonequilibrium systems. Indeed if a macroscopic system 
is maintained out of equilibrium by time independent forces it is not uncommon that it 
exhibits oscillations, which may be periodic or nonperiodic (chemical oscillations, hydro- 
dynamic turbulence). A time independent state p describing such an oscillating system is 
expected to have a natural decomposition 



where (/*)*Pa = Pe l \ and p is ergodic under (#*). It is however not clear how to express 
such a decomposition at the microscopic level. 

3. Thermodynamic limit. 

Statistical mechanics deals mostly with macroscopic systems, and some important 
phenomena like phase transitions make sense only in the large system limit. At a technical 
level, an important feature of the large system limit (thermodynamic limit) in equilibrium 
statistical mechanics is the equivalence of ensembles. Roughly speaking, this corresponds 

* One can make U 1-valued by going to the covering space of phase space, this approach 
is taken by Hoover et al. [52], Dettmann and Morriss [23]. 

** As an example, if the Anosov flow (/*) is not mixing, i.e., is a suspension of an Anosov 
diffeomorphism /, one can write 



= lim (f)*m 



(1) 



t— >oo 





where p is an /-ergodic measure on a manifold M with dimM = dimM — 1. Here (f T )*p 
depends periodically on r. 
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to the fact that for large systems, extensive variables usually fluctuate very little. One 
example, which will be of interest to us, is that if one fixes the (total) energy, then the 
(total) kinetic energy fluctuates very little (with respect to the microcanonical ensemble, 
i.e., the normalized Liouville measure on the energy shell). This is really a fact about 
geometry in a large number of dimensions** which is reasonably well understood in the 
situation of equilibrium statistical mechanics, thanks to the convexity properties of the 
entropy (see Ruelle [78], Lanford [59]). The notion that extensive variables often do not 
fluctuate much in nonequilibrium statistical mechanics (with the microcanonical ensemble 
replaced by the SRB measure) probably has some validity, but a firm theoretical basis is 
lacking. For a physical discussion see the papers of Evans and Sarman [33]*, Cohen and 
Rondoni [19]. 

There are also some danger signs coming from physics. Suppose for instance that we 
try to produce a nonequilibrium steady state corresponding to homogeneous shear in a 
fluid. Macroscopic fluid theory tells us that the homogeneous shear flow that we try to 
produce is unstable and will become turbulent with, presumably, fluctuations on all scales. 

4. Dimensionality loss and the definition of entropy. 

Let Ai > A 2 > • • • > Xd be the Lyapunov exponents associated with the SRB measure 
p. If d is the largest integer such that 



the Lyapunov dimension is defined to be 



X 



d+l | 

by Kaplan and Yorke [54], [35]. According to these authors, Dl is in many cases equal to 
the Hausdorff dimension of p. 

For a macroscopic system outside of equilibrium it is expected that the Lyapunov 
exponents have an asymptotic density 0(-) when D (the total dimension of phase space) 
tends to infinity. The number of Lyapunov exponents in the interval (Ai,A2) is thus 
D f* 2 (j)(X)d\. It is usually assumed that the support of 4>(-) is an interval [A m ; n , A max ] 
where A m ; n , A max remain finite for D — > oo (in fact A m i n , A max may diverge exponentially, 
see Searles et al. [89]). The entropy production is 

e = —D / <j)(X)XdX (2) 



** A typical fact of this sort is that in many dimensions, most of the mass (=Lebesgue 
volume) of a solid sphere is concentrated near its surface. 

* This paper gives a perturbative treatment of the equivalence of the isokinetic and 
isoenergetic constraints, but ignores the singular character of the measures to be used. 
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For a Hamiltonian system we have 0(A) = 0(— A), while if pairing holds we have 0(— A — 
a) = 0(A) where Da/2 = e. 

To obtain the Lyapunov dimension Dl of p, we determine Ao such that 

P A ma x 

/ 0(A)AdA = O (3) 



and write 

The dimension loss is thus 
with 



D L =D 0(A) dA 



A = D - D T . = DS 



Ao 

0(A) dX (4) 



Note that A is an extensive variable {i.e., proportional to the number of degrees of freedom 
of the system). From (2) and (3) we obtain 



e f x ° 
~P = / 0(A) A 

J A m i n 



Suppose that A m ; n and A max remain finite when D — > oo. Close to equilibrium we have 
A ~ A min w -A max hence (4) gives 



----A 

j-y A max 

and therefore 

A = D5 



Amax 

(see Posch and Hoover [73]). 

Corresponding to a nonequilibrium steady state p with A > we have a zero D- 
volume of occupied phase space, and therefore an entropy — oo. We may however try to 
define a "Lyapunov entropy" corresponding to -D^-volume in phase space. Suppose thus 
that V(e, d) is the occupied volume in phase space* in dimension d associated with the 
length scale e. We refrain from trying to give a precise definition of V(e, d), but we assume 
a scaling form 

V{e, d) w Ce d ~ DL 

for small e (such a form is familiar, for instance from renormalization group theory). 
Therefore the specific entropy 

— \ogV{e, d) w — logC+ — - — loge 



* In the case of n identical particles we include as usual a factor 1/n! 

60 



depends weakly on e if d is close to D L , and the "specific" Lyapunov entropy D 1 logC is 
reasonably well defined if one takes e 2 of the order of Planck's constant. 

One way to express the above formulae is that in our nonequilibrium situation a 
fraction 6 of the degrees of freedom is frozen, and that the Lyapunov entropy Sl = logC 
is the entropy corresponding to the fraction 1 — 6 that is not frozen. But can one give a 
physical meaning to Sx? 

We have up to now considered classical systems. Let us now assume for the sake of the 
discussion that our description is a classical approximation applied to a quantum system. 
According to this bold assumption, the "Lyapunov entropy" Sl is now the logarithm of a 
number of available degrees of freedom*, and is thus in principle a physically meaningful 
quantity. (Note that in the quantum accounting, the frozen degrees of freedom make an 
additive contribution equal to 0, which can be ignored, while in the classical situation, the 
additive constant is — oo, which cannot be forgotten). 



* The quantum description fixes a unit of volume in phase space, based on Planck's 
constant, and eliminates the usual arbitrary additive constant in the definition of the 
entropy. 
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